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New thermoelectric materials with superior transport properties at high temperatures have been discovered. 
These materials are part of the large family of skutterudites, a class of compounds which have shown a 
good potential for thermoelectric applications. The composition of these novel materials, called filled 
skutterudites, can be represented by the formula LnT,Pn,, (Ln = rare earth, Th; T = Fe, Ru, Os, Co, Rh, Ir; 
Pn = P, As, Sb). In these compounds, the empty octants of the skutterudite structure which are formed in 
the TPn, (~ T,Pn,,) framework are filled with a rare earth element. Some of these compositions, based on 
CeFe,Sb,,, have been prepared and have shown exceptional thermoelectric properties in the 350-700°C 
temperature range. At room temperature, CeFe,Sb,, behaves as a p-type semimetal, but with a low thermal 
conductivity and surprisingly large Seebeck coefficient. These results are consistent with some recent band 
structure calculations on these compounds. Replacing Fe with Co in CeFe,Sb,, and increasing the Co:Fe 
atomic ratio resulted in an increase in the Seebeck coefficient values. Measurements on bulk samples with 


a CeFe, ;Co, ;Sb,, atomic composition and p-type conductivity resulted in dimensionless figure of merit 
ZT values of 1.4 at 600°C. 
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Thermoelectric Science 
and Technology - An 


Overview 


John C. Pazik, Guest Editor 
Office of Naval Research 


Thermoelectric (TE) materials, and devices prepared 
from them, are of interest to the Navy for two basic reasons: 
1) the devices can be used in reliable, quiet, solid state cool- 
ers and 2) the devices can be used for exhaust heat recovery 
and/or power conversion. Based on these two fundamental 
capabilities (cooling and power generation), thermoelectric 
materials can have broad impact on Navy warfighting capa- 
bilities including Forward Presence, Readiness, Support and 
Infrastructure, and Joint Strike. Exhaust heat recovery and 
power generation aspecis of thermoelectric materials may 
ultimately lead to reduced fuel usage which “can minimize 
the logistic footprint without loss of capability” and “increase 
naval facilities’ energy independence.” Since TE devices 
used for cooling applications are solid state (no moving parts, 
lubricants, seals, etc) they also impact the area of “design-in 
reliability, maintainability, and availability.” From an envi- 
ronmental perspective, the non-CFC based cooling offered 
by TE technology can positively impact the Navy’s “capa- 
bility to operate unimpeded by environmental issues” and 
provides a “substitute for hazardous materials and processes.” 
Ultimately one can envision TE cooling devices which “al- 
low strike assets to manage their signatures to delay detec- 
tion and targeting from threat sensors.” Specific cooling 
applications where TE technology is currently used in Navy 
systems (or envisioned for use) include ship-board air-con- 
ditioning, IR detection systems, cryogenic electronics (mi- 
croprocessors, high temperature superconductors), and en- 
vironmental suits. Thermoelectric power generation appli- 
cations include exhaust heat recovery (gas-turbine engines, 
exhaust stacks, diesel engines) and silent auxiliary power 
sources (submarines, UUV). Commercial applications of 
TE technology are equally enormous and include cooling 
units for cellular base stations, high technology home refrig- 
eration, increased fuel efficiency for trucks and automobiles, 
and recreational coolers. 

Clearly there are other technologies available which 
can meet many of the cooling and power conversion needs 
of the Navy and private sector. So what exactly is it about 
thermoelectric technology that makes it so attractive for the 
various applications? The simple answer is that TE cooling 
and power generation devices are solid state. There are no 
moving parts, no refrigerants, no lubricants, and no seals. 
Quiet, non-emissive, vibrationless, and reliable operation are 
the key aspects of TE technology which set it apart from all 
the competing technologies. If TE technology is so great 
why isn’t iteverywhere? Quite simply, the efficiency of cur- 
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rent systems is too low to be used in anything but the most 
demanding applications where the higher cost (i.e. low effi- 
ciency) is outweighed by the mission requirements. For ex- 
ample, TE has been the technology of choice for power gen- 
eration on NASA deep space missions because of its inher- 
ent long term reliability and vibrationless operation. 

The most critical issue for researchers in TE science 
and technology today is the development of new materials 
and devices with intrinsic properties much better than those 
currently in use today. It’s helpful at this point to recall the 
basic thermoelectric effects. When any material is placed in 
a temperature gradient a potential will be produced which 
will depend on the magnitude of the temperature gradient 
and the physical properties of the material. Similarly, if a 
current is passed through the junction of two dissimilar ma- 
terials, heat will either be liberated or absorbed at the junc- 
tion depending on the direction of current flow and material 
properties. The magnitude of these “thermoelectric effects” 
determine the quality (figure of merit) of a given set of ther- 
moelectric materials. A typical thermoelectric device on the 
market today is based on semiconductor materials which are 
n and p doped and connected in series (as depicted on the 
front cover of this issue). For a cooling application, one can 
conceptually imagine that electrons will gain or loose en- 
ergy (depending on direction of current flow) to the surround- 
ings as they move over the potential barrier at the p-n junc- 
tion. The ultimate performance of a thermoelectric based 
system (cooling or power generation) is determined by the 
intrinsic materials figure of merit (FOM) defined in equa- 
tion | where S is the Seebeck coefficient, 6 is the electrical 
conductivity, and x is the thermal conductivity, and T is the 
absolute temperature. 


S°oT 
K 


ZT = (1) 
Equation | provides general guidance for the type of materi- 
als needed to maximize ZT and ultimately improve device 
performance. The Seebeck coefficient is a measure of the 
ability of a material to generate an EMF in a temperature 
gradient. The larger the EMF the material can generate for a 
given temperature gradient the better the material will be for 
TE technology. It is also desirable to maximize the electri- 
cal conductivity, 6, because a large 6 will result in small 
internal heat loss (PR heating). Finally, one wishes to choose 
a material which has low thermal conductivity, K, since this 
will minimize the equalization of the temperature gradient. 
One quickly sees that the relationship between the various 
material parameters is complex and must be considered when 
designing new materials. For example, « has two compo- 
nents, an electrical and phonon contribution (equation 2) and 
the electrical component, K_,, is related to oO via the 
Wiedemann-Franz Law. 


el’ 


K=K,+K,, (2) 





Essentially any increase in 6 will result in an increase in kK 
(via ,,) and thus the effective ZT will not be improved. Thus 
for a material whose electrical component is the major con- 
tribution to the thermal conductivity (i.e. a metal), strategies 
to improve ZT via improvements in electrical conductivity 
or reduction in thermal conductivity will not be fruitful. 
Metals typically have low Seebeck coefficients (1-10 
ut'V/K) and high electrical (210° ohm-cm") and thermal con- 
ductivity (10-1000 W/m-K) giving a relatively low FOM, 
~10°. Conversely, insulators have high Seebeck coefficients 
(1000 uL'V/K) and low thermal conductivities (.01 -1 W/m-K), 
however the electrical conductivity is extremely low 
(<10* ohm-cm") resulting in ZT’s around 10°. Semicon- 
ductors are currently the materials of choice. Seebeck coef- 
ficients range from 100-500 y'V/K, electrical conductivity 
ranges from 10°-10* ohm-cm", and thermal conductivities 
between 1-100 W/m-K. The best, well understood, thermo- 
electric materials' are based on the Bi-Te-Se-Sb alloys and 
have FOM’s of approximately 1.0. It is also important to 
recognize that the fundamental materials properties are tem- 
perature dependent. So a material may make an outstanding 
thermoelectric device over one range of temperatures but 
could and is likely to be useless over other temperature ranges. 
Si-Ge materials function optimally between 1100-1300K 
while Bi,Te, based materials operate best between 200-500K. 
A final word about measuring the performance of ther- 
moelectric materials and devices. In addition to the intrinsic 
FOM, the TE device performance properties may be defined 
by a coefficient of performance (COP), equation 3 and 4 


Cooling Equation: 


(fivaT )-— 
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Power Generation Equation: 
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where T. is the temperature of the cold side and T,, is the 
temperature of the hot side. It is important to realize that 
ultimately one must consider the AT available for power gen- 
eration or for the cooling requirement and that the actual 
device efficiency is dependent on this temperature differen- 
tial. When AT is small, TE based cooling technologies be- 
come more competitive than traditional Rankine cycle sys- 
tems for a given cooling application. It should also be pointed 


out that for ZT values ranging from 3-8, TE based cooling 
could compete head on with Rankine cycle cooling in al- 
most every cooling application. 

This issue of Naval Research Reviews provides an 
overview of the direction current thermoelectric research is 
taking. The five general articles represent the core of the 
ONR program in thermoelectric materials science and tech- 
nology: theory, experiment and property measurement. The 
article by Singh and Reinecke (et. al.) describes the current 
theoretical approaches to designing better thermoelectric 
materials. It includes discussions of band structure calcula- 
tions for complex inorganic materials and describes meth- 
ods of improving ZT through quantum confinement. 
Venkatasubramanian has provided a report on experimental 
efforts currently underway to synthesize thin film TE struc- 
tures and how the design of these structures could impact 
the FOM. Kanatzidis and DiSalvo provide an overview of 
current approaches in solid state chemistry for synthesizing 
bulk materials which may exhibit enhanced ZT. The focus 
of their discussion involves molten flux synthesis of multinary 
chalcogenides and an exploration of valence fluctuation com- 
pounds, metals with unusually large Seebeck coefficients. 
Slack and Fleurial describe an interesting class of compounds 
referred to as skutterudites'. These materials have large, 
complex unit cells and have properties which show promise 
as TE materials. Uher provides an overview of techniques 
currently being used to measure the properties of new ther- 
moelectric materials. This is a critical area of research since 
accurate, reproducible measurements are critical to the iden- 
tification of promising materials. 


‘Recent reports indicate that a new material of nominal composition, CeFe, Co, Sb,,(0 < x < 1), has a FOM of 1.4 (see inside front 


cover). 
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Theoretical 


Approaches for 
Understanding and 
Optimization of 
Thermoelectric 


Materials 


D.J. Singh, T.L. Reinecke, W.E. Pickett and P.J. Lin-Chung, Naval Research Laboratory, Washington, D.C. 


Introduction 


Materials research has unearthed a plethora of inter- 
esting, useful and often unanticipated phenomena. Some 
more notable recent examples include the fullerenes, high 
temperature superconductivity, colossal magnetoresistance 
in manganites, the quantum hall effects, superhard magnetic 
materials, ductile ordered structural intermetallics and oscil- 
latory couplings in magnetic multilayers. Discoveries have 
often come through serendipity or at least through a great 
deal of searching based on physical intuition and qualitative 
theoretical ideas. No doubt there are many materials with 
unanticipated but useful properties remaining to be found. 
However, the task of finding a material meeting a given set 
of specifications is often akin to finding the proverbial needle 
in a haystack. This article describes advances in theory and 
modeling that permit more systematic approaches and illus- 
trates these as they are being applied to ONR supported re- 
search seeking better thermoelectrics (TE). Two aspects are 
highlighted: (1) The use of first principles calculations to 
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understand materials properties and thereby point to possible 
improvements; and (2) The use of theory to understand 
nanostructural effects and to suggest novel quantum well 
structures for enhancing TE properties. These comprise the 
two main sections of this article. 

Recent years have seen tremendous advances in our 
ability to fabricate and characterize a wide variety of tailor- 
made structures on nanometer (nm) length scales, using such 
techniques as molecular beam epitaxy, organo-metallic va- 
por phase epitaxy, lithography and etching as well as novel 
chemical processes. The properties of systems structured on 
these small length scales are often radically different from 
the bulk materials of which they are composed and depend 
on a wide range of changeable variables such as specific 
materials compositions and well widths. 

Fabricating a sufficient variety of materials and struc- 
tures using the Edisonian approach for testing can be time 
consuming and costly. It has been found in this field that 





model calculations of the properties of these systems carried 
out in conjunction with experiment and fabrication are a 
powerful tool in predicting novel phenomena, in understand- 
ing observations and in guiding the fabrication of new mate- 
rials and structures. The electronic states involved in these 
systems are often shallow states near band edges, and the 
length scales of interest typically are of the order of 10’s of 
nanometers or more. Under these conditions the electronic 
properties can be described very well within generalized ef- 
fective mass methods and transport within Boltzmann equa- 
tion approaches. 

Searching for improved materials by synthesizing and 
characterizing great numbers of different compositions and 
structures, guided by physical intuition and simple models, 
has lead to remarkable discoveries. However, recent advances 
in theory, both of the computational variety driven by new 
algorithms and data processing technology, as well in model 
calculations for structures of experimental interest, have pro- 
vided alternate and more efficient strategies. These are itera- 
tive in nature and in the best developed cases consist of 
strongly interacting theoretical, computational and experi- 
mental efforts. Theory, both formal and computational, play 
several roles in these efforts. Besides suggesting promising 
candidate materials and structures, theory plays a key role in 
unraveling the properties of new materials, and the under- 
standing thus gained provides important guidance in the sub- 
sequent optimization of the properties. 

The work described here is aimed at developing im- 
proved thermoelectrics. TE materials are characterized by 
the TE figure of merit Z = oS?/k, where o is the electrical 
conductivity, S is the Seebeck coefficient (also called 
thermopower) and x is the thermal conductivity. ZT (T is 
temperature) is a dimensionless quantity which slightly ex- 
ceeds one at present in the best materials. Higher values of Z 
yield better TE performance. There is no known fundamen- 
tal limit on Z; Carnot efficiency corresponds to infinite Z. 


Understanding Materials 
from First Principles 


Materials are fundamentally collections of electrons 
and atomic nuclei; the interactions are those of charged par- 
ticles (Coulombic) and the dynamics are those of quantum 
mechanics. From this reductionist perspective, the underly- 
ing principles have been well understood since the early part 
of the 20th century. Nonetheless, most important discover- 
ies in materials science were unanticipated and several have 
proven remarkably resistant to explanation from basic prin- 
ciples. Clearly, this perspective misses the tremendous com- 
plexity and richness of phenomena that is materials science. 
Substantial progress in understanding materials from first 
principles has only occurred in the last decade or so. The 
origin of the difficulties in bridging the microscopic physics 


and the real world lies in the quantum mechanical equations 
and the large number of particles. Although simple in form, 
these are 3N dimensional differential equations, where N is 
the (practically infinite in materials) number of particles; as 
a result the quantum mechanics is intractable in its most ba- 
sic form. 

The major breakthrough that allowed progress from 
first principles came in the mid 1960’s from ONR sponsored 
basic research in quantum mechanics, and consisted of the 
development of density functional theory (DFT) by 
Hohenberg, Kohn and Sham [1] and especially practical com- 
putational schemes based on it. The DFT and particularly 
the local density approximation (LDA) to it, along with al- 
gorithmic developments and advances in computing tech- 
nologies, have now made possible realistic studies of a wide 
variety of materials and properties. The DFT framework re- 
duces the intractable problem of solving of a 3N dimensional 
differential equation to that of solving a set of N three di- 
mensional equations, albeit coupled by a non-linear, density 
dependent potential. With modern computers and algorithms, 
representative samplings of these systems of equations can 
be solved on a routine basis, provided that the complexity of 
the material, as characterized by the unit cell size, is not too 
great. 


Computational Approaches for 
Thermoelectrics 


The performance of a thermoelectric material is char- 
acterized by the figure of merit, Z. The numerator of Z con- 
sists of electronic transport quantities; these depend on the 
band structure (the energies of allowed states as functions of 
their net crystal momentum) near the Fermi energy and the 
scattering mechanisms for the charge carriers. These terms 
depend strongly on the doping level with the result that ob- 
taining high Z requires careful optimization of the doping 
level, and in fact even discerning whether or not a given semi- 
conductor has TE potential requires exploration of the varia- 
tion with doping. Knowledge of the band structure, which 
can be obtained reliably using modern electronic structure 
methods, permits the doping dependence to be calculated, 
thereby bypassing time consuming experimental searches. 

The denominator, «, is the sum of an electronic ther- 
mal conductivity (related to 6) and a lattice contribution, the 
latter normally dominating in semiconductors. Understand- 
ing K requires an understanding of the lattice vibrations and 
the scattering processes affecting them. 

The work discussed here is based on the Linearized 
Augmented Planewave Method [2], which is very general 
with regard to the materials that it can be applied to and is 
state-of-the-art as regards accuracy in solving the DFT equa- 
tions. No adjustment of the DFT band gaps is made. DFT 
gaps, although smaller than experiment in certain semicon- 
ductor classes, agree well with experiment in well hybrid- 
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ized transition metal compounds, like skutterudites. Because 
these methods are essentially microscopic and derived from 
the underlying quantum mechanics and interactions, they do 
not depend on experimentally determined parameters. Rather, 
one may calculate properties of a material, starting from just 
the atomic numbers and usually the crystal structure. This 
first principles aspect is often crucial in understanding novel 
materials, as experiments are often unavailable, incomplete 
or subject to interpretation. The importance of such an ap- 
proach is well illustrated by the electronic structures of the 
binary skutterudites, which have been under investigation as 
potential TE materials. Theory has led to a complete revi- 
sion in the description of their TE properties. 


Binary Skutterudites: New 
Excitement from Old Materials 


The binary skutterudites, MX, (M = Co, Rh, Ir; X = P, 
As, Sb), although cubic in symmetry, have a rather complex 
crystal structure (Fig. 1) containing sixteen atoms per cell 
and large open voids. They have been known and generally 
characterized as rather conventional semiconductors with 
moderate band gaps [3], although Ackerman and Wold did 
observe a weak optical signal extending to low energy and 
questioned this characterization. Recently, these materials 


have come under investigation as a potential basis for im- 
proved thermoelectrics [4,5]. 

Slack and co-workers measured high carrier mobili- 
ties (i.e. high o for a given doping level) and reasonably large 
values of the Seebeck coefficient in hole doped IrSb,. How- 
ever, the measured kK was more than an order of magnitude 
higher than the estimated minimum thermal conductivity (see 
the article by Slack in this issue), resulting in relatively low 
values of Z at room temperature. Furthermore, since TE is a 
bulk application, the use of Ir (an expensive precious metal) 
raises serious cost issues. Several questions were raised by 
these seminal studies: (1) What is the band structure and what 
does it imply for the variation of Z with doping level and 
temperature; specifically, is there a doping regime where this 
material has a high ZT? (2) What modifications can be made 
to reduce « while maintaining electrical properties favorable 
for high ZT? (3) Can Ir be replaced by any less expensive 
elements without degrading TE performance? 

The starting point for addressing these questions is the 
electronic band structure, which we calculated from first prin- 
ciples [6]. Accurate calculations for complex materials such 
as these were made possible by use of the supercomputers 
supported by the DoD high performance computing mod- 
ernization program and algorithmic improvements in the 
LAPW method made with ONR support at NRL and else- 
where. 





Figure 1. 


Crystal structure of the binary skutterudites, MX,. This body 
centered cubic structure may be viewed as consisting of 
interpenetrating simple cubic lattices, the first consisting of 
M atoms on each lattice site and the second consisting nearly 
square X, rings oriented with alternating x, y, z and empty 
sites. The filled skutterudite structure derives from filling these 
empty sites with a third element. First principles methods 
described in this article allow computation of electronic, 
vibrational and other properties from knowledge of the crystal 
structure. 
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Figure 2. 


Electronic band structure of IrSb, along representative 
direction in the Brillouin zone. The dashed horizontal line 
denotes the Fermi energy (E,). Note the valence band 
dispersing upwards to the conduction band minimum at the 
I point. This is at odds with previous views of the material as 
a moderate band gap semiconductor. 








The band structure of IrSb, is shown in Fig. 2, and the 
corresponding electronic density of states (DOS) in Fig. 3. 
The apparent complexity of these plots results from the fact 








Figure 3. 


Electronic DOS of IrSb, and projections. E, is denoted by the 
dashed vertical line. Note the low DOS in the gap-like region 
just below E.. 
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that the number of bands is roughly proportional to the num- 
ber of atoms in the unit cell (2 for Silicon or GaAs; 16 for 
IrSb,). There is a rather well defined pseudo-gap between 
manifolds of valence (occupied states below the Fermi en- 
ergy, E,) and conduction (unoccupied states above E,) bands, 
and this is reflected in the DOS by very low values in the 
first 1.2 eV below E,. This corresponds to the gap seen in 
optical experiments and assigned as the band gap. 
However, the DOS in this region, although small, is 
not exactly zero, and this non zero DOS is critical to under- 
standing the electrical transport (6 and S). Instead, there is a 
single band that crosses this pseudogap reaching the con- 
duction band minimum at the [ point. The reason why this 
was not previously established in experiments is that it is 
one dispersive band out of many, and therefore gives a very 
weak optical signal that is difficult to distinguish from arti- 
facts e.g. due to defect states in imperfect samples. All re- 
ported measurements on IrSb, were on lightly hole doped 
material, which apparently forms most readily. In this case, 
E, lies slightly below the nominal (undoped) position, and 
the electrical transport properties are determined almost en- 


tirely by the properties of this unanticipated band. 

In ordinary semiconductors, the valence and conduc- 
tion bands are parabolic in shape at least near E,, so that the 
band is a smooth curve reaching (for valence bands) its maxi- 
mum with zero slope. The transport properties are then un- 
derstood in terms of the effective mass m*, which gives the 
curvature, and a handful of other parameters, such as scat- 
tering times. Remarkably, the valence band in IrSb, does not 
fit this picture. Rather than being parabolic as in conven- 
tional band structures, it is linear almost exactly up to the 
band edge. For hole doping above about n = 3x10"* cm? (TE 
materials typically have n ~ 10"° cm*), E, will occur in this 
linear regime with the implication that the electrical trans- 
port will scale differently than in conventional semiconduc- 
tors. For example, S is predicted to vary as n™? instead of 
n*? and o will vary approximately as n”* instead of n. Our 
calculations show a very similar linearly dispersing band 
crossing the gap in CoSb,; CoAs, also has a gap crossing 
state but its dispersion has the conventional parabolic shape. 
(CoAs, becomes interesting because it is predicted to be a 
member of the very small class of true gap-less semiconduc- 
tors.) These predictions, which require a revision of the pre- 
vious understanding of the materials, are in quantitative agree- 
ment with the Seebeck coefficients measured by Slack and 
have now been confirmed by additional experiments testing 
these scalings and carefully searching the gap region [7,8]. 

Besides the new semiconductor physics that this dis- 
covery opens up, there are important implications for TE: 
(1) Zin p-type IrSb, will depend only weakly on doping, so 
high Z will not arise in p-type IrSb, through optimization of 
n alone; (2) the linear dispersion yields large shifts in E, for 
given n, so non-degenerate behavior and roughly linearly 
increasing S will persist to high T (favorable for high tem- 
perature TE applications); (3) the parabolic conduction bands 
have reasonably high effective masses (~ 0.5 m,), so that if 
n-type mobilities are comparable to p-type values, the n-type 
material will have significantly higher ZT after optimization 
of n; and (4) highly covalent electronic band structures are 
found including notably strong bonding in the Sb, rings. The 
implication of this last point is that « in CoSb, should not be 
much degraded relative to IrSb,. Since the electronic proper- 
ties are similar, the cheaper Co based materials may be good 
substitutes for Ir containing skutterudites. 


Evaluation before Synthesis 


These studies of the binary skutterudites indicated that 
modifications that would lower « and/or improve the electri- 
cal properties are needed to make them into competitive room 
temperature TE materials. This suggests studies of n-type 
binary skutterudites, which may have improved electrical 
properties (numerator of Z), and the exploration of strate- 
gies for lowering the thermal conductivity (denominator of 
Z). 
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Slack and co-workers noted that « of IrSb, is more 
than an order of magnitude higher than theoretical estimates 
of the minimum « for skutterudite and that the structure con- 
tains large voids that should accept atoms. Based on this, he 
speculated that x could be lowered considerably if weakly 
bound atoms could be inserted into these voids, where they 
might “rattle” about and scatter low energy phonons. He sug- 
gested the possibility of inserting heavy rare gas atoms, par- 
ticularly Xe, by high pressure synthesis. This choice was 
motivated by the expectation that these inert atoms would 
not strongly scatter electrons, presumably leaving the elec- 
tronic properties unchanged, and that being weakly bound 
they might be particularly good scatterers of low energy 
phonons. However, testing these conjectures experimentally 
is difficult because of the need to construct a high pressure 
apparatus, the high cost of Xe gas, and the potential need to 
explore a wide range of processing conditions. Accordingly, 
we used first principles calculations for the postulated com- 
pound, Xelr,Sb,, to explore some of these issues. 

The calculated band structure of Xelr,Sb,, was found 
to be practically identical to that of the corresponding un- 
filled skutterudite, IrSb,, with no discernable change near 
E,. This lends strong support to one of Slack’s conjectures, 
i.e. Xe will not scatter electrons so that the electronic prop- 
erties will not be degraded relative to IrSb,. Unfortunately, 
this is not the whole story. 

The DFT computational framework directly gives the 
total binding energies of solids. These can be compared to 
obtain the relative stability of various structures, as well as 
to compute such properties as defect energies, diffusion bar- 
riers and phonon spectra. Calculations for Xelr,Sb,, indi- 
cate that the compound is energetically unfavorable relative 
to separated Xe and IrSb, by 5-10 Kcal/mole. This substan- 
tial energy, and the inertness of Xe, which makes chemical 
synthesis routes unlikely, suggests that Xelr,Sb,, will be ex- 
ceedingly difficult to stabilize, if it can be done at all. Inser- 
tion of lighter (and smaller) rare gas atoms is more likely to 
be feasible. 


Alternate Filled Skutterudites 


Although it may be difficult to synthesize rare-gas filled 
skutterudites, the basic idea of filling the voids in the struc- 
ture to obtain lower thermal conductivity may be sound, and 
in fact a number of filled skutterudite structure compounds 
are known to exist. These are derived from the binary 
skutterudites (e.g. CoSb,) by replacing the transition metal 
by the element to the left in the periodic table, and inserting 
a di- or tri-valent cation forming e.g. BaFe ,Sb,, or LaFe,Sb... 
These compounds have lower valence electron counts than 
the binary skutterudites since replacement of the transition 
metal results in four fewer electrons (one for each of four 
atoms), and the cation insertion leads to the addition of two 
or three valence electrons. Compounds without the transi- 
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tion metal substitution would have higher electron counts, 
but these have not been successfully synthesized, a fact prob- 
ably related to the tendency of binary skutterudites to form 
with p-type (lower electron count) dopings and the 
antibonding nature of the conduction band states. 

Metallic conduction is generally unfavorable for TE 
applications because normal metals have Fermi energies ly- 
ing in broad bands leading to very low Seebeck coefficients. 
(Note the exception of f-band valence fluctuation materials 
discussed in this volume by DiSalvo.) Unfortunately, almost 
all of these filled skutterudites are known to be metallic. 
LaFe,P,, is even a superconductor. 

This metallic character has been rationalized by sup- 
posing that the band structures are substantially the same as 
in the binary skutterudites, but that the Fermi level is low- 
ered corresponding to the reduced electron count. Within this 
so-called rigid band picture, dispersive bands will cross the 
Fermi level and the material will be metallic. However, there 
is no good reason for supposing that these materials follow 
rigid band behavior, and in fact there are some semiconduct- 
ing filled skutterudites, notably UFe,P,,, CeFe,As,, and 
CeFe,P,,. CeFe,Sb,, has been reported as a possibly low car- 
rier density metal, although it is unclear whether this is due 
to experimental difficulties in finding a very small gap in 
doped samples. 

Semiconductivity in CeFe,X,, (X=P,As) has been ra- 
tionalized in terms of the rigid band picture (above) by sup- 
posing that Ce is tetra-valent ia these compounds, yielding 
the same electron count as CoX,. However, the fact that the 
much larger trivalent La ion enters the same structure, with 
the implication that tri-valent Ce would fit as well, is at odds 
with this view. 

We have used first principles calculations of the elec- 
tronic structure to resolve this issue and to assess the TE 
potential of the Ce filled skutterudites. One immediate diffi- 
culty was that the crystal structures (normally an input to the 
calculations) of CeFe,As,, and CeFe,Sb,, were only partially 
known, and in particular the two parameters giving the 
pnictogen (P, As, Sb) positions in the unit cells had not been 
determined experimentally. At the same time, Morelli and 
Meisner [9] reported transport measurements on CeFe,Sb.,, 
raising questions about its electronic structure and showing 
that it does indeed have a much lower thermal conductivity 
than CoSb, - a result that increases interest for potential TE 
application. As mentioned, one of the main quantities that 
density functional calculations provide is the total energy 
for any given structure. Accordingly, it is possible to deter- 
mine structural parameters by performing calculations of the 
energy as a function of these coordinates, thereby determin- 
ing the minimum energy configuration, which then corre- 
sponds to the actual crystal structure. A by-product of this 
mapping of the energy surface is its curvature about the mini- 
mum, which yields vibrational frequencies or phonons. This 
prescription was used to determine the pnictogen positions 
and the associated A, Raman active phonon frequencies in 





CeFe,As,, and CeFe,Sb,,. A full set of phonon dispersions 
was calculated using this procedure combined with force 
constant models for the unfilled skutterudite CoSb, [10]. 

Our first principles calculations for CeFe,P.., 
CeFe,As,, and CeFe,Sb,, show that these materials have 
qualitatively similar electronic structures with smooth varia- 
tion from the lighter to the heavier pnictogen compounds 
[11]. All three compounds are found to be semiconducting, 
albeit with small gaps of 0.34, 0.16 and 0.10 eV in the P, As 
and Sb compounds, respectively. The band structure for 
CeFe,As,, is shown in Fig. 4. 

The topmost valence (occupied) bands consist of hy- 
bridized Fe 3d and pnictogen p states, while the lower va- 
lence bands are mainly dominated by pnictogen p character. 
On the other hand, the lowest conduction bands are com- 
prised of nearly flat spin orbit split manifolds of Ce 4f states, 
and are apparently situated in a gap formed by the other bands. 
However, analysis of the wavefunctions reveals that this pic- 
ture is only qualitative. 

In fact, there is considerable Ce 4f character in bands 
below E,. When integrated, this occupied 4f character 
amounts to one electron per Ce atom, i.e. tri-valent Ce at 
odds with the previous explanation of the semiconducting 
character as being due to tetra-valency of Ce. Rather, the 
band gaps in these compounds arise because of hybridiza- 
tion or mixing of Ce 4f states with the skutterudite Fe 3d - 
pnictogen p states. This hybridization, which is strongest in 
the phosphide and weakest in the antimonide, leads to an 
interesting combination of small band gaps (possibly over- 
estimated for f-bands within the LDA), and large effective 
masses. This is in contrast to the low effective masses more 
commonly found in small band gap materials. 





Figure 4. 


Band structure of CeFe As,.. The dashed horizontal line 
denotes E.. 








High band masses lead to high Seebeck coefficients 
for a given doping level, and as such are favorable for TE 
applications, provided other conditions are met as well. In 
particular, a low lattice thermal conductivity, K, is needed, 
along with reasonably high carrier mobilities (for high ©). 
High mobilities arise in materials with highly itinerant band 
states. In these materials, the hybridization, and itinerancy is 
strongest in CeFe,P,,, and weakest in CeFe,Sb,,. The ars- 
enide is intermediate. 

Morelli and Meisner have, as mentioned, shown that 
CeFe,Sb,, does indeed have very low values of «. Calcu- 
lated Raman phonon modes in CeFe,As, , are similar to those 
of CeFe,Sb,,, but are shifted to about 50% higher frequency. 
Based on this and the highly covalent picture of the bonding 
that emerges from analysis of the calculated electronic struc- 
tures, it seems likely that CeFe,As,, also has low x, though 
probably not as low as the antimonide. The presence of light 
atoms in the phosphide and its considerably contracted lat- 
tice parameter, which would not favor soft Ce modes, imply 
that its thermal conductivity may be substantially higher. 
These considerations, and the calculated band structures, 
suggest that the arsenide and antimonide may be favorable 
TE materials. Experimental measurements of the transport 
properties of n- and p-type crystals will be needed to estab- 
lish this. 


Nanoscale Structures: 
High Performance 
Superlattices 


Remarkable advances have been achieved in recent 
years in our ability to fabricate a wide variety of novel mate- 
rials in the form of superlattices and other low dimensional 
and composite structures. Among the best known of these 
structures are quantum well superlattices, prepared in the 
form of multilayers. The existence of conduction and va- 
lence band offsets between adjacent material layers creates 
potential wells separated by potential barrier layers. The 
charge carriers in the quantum wells are confined to move 
perpendicular to the superlattice direction and thus give rise 
to two-dimensional conduction. Recent advances in growth 
and microfabrication techniques, as well as etching and li- 
thography techniques, have also made it possible to obtain 
quasi-one-dimensional wire like superlattice structures of one 
material within another. Such wires have cross-sectional ar- 
eas ranging from a few atoms to square micrometers or more. 

The properties of these systems can be modified in 
controlled ways by a variety of changes in structure and 
materials. For example, the potential off-sets between dif- 
ferent materials can give rise to the confinement of carriers 
in effectively two dimensional or in one dimensional sub- 
bands. In two dimensions the density of states is proportional 
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to the inverse of the well width, and in one dimension to the 
inverse square of the effective wire width. This can lead to 
greatly enhanced densities of states and transport properties. 
In addition, the thermal properties of these systems can be 
changed through changes in the phonon scattering at inter- 
faces and due to disorder (Recall that decreases in the ther- 
mal conductivity lead to more desirable TE materials). For 
reasons such as these, superlattices and related composite 
systems are of considerable interest in the search for materi- 
als with improved thermoelectric properties. 

An interesting proposal involving superlattice systems 
for TE applications was made recently by Hicks and 
Dresselhaus [12] who argued that the TE figure of merit of a 
material, for transport parallel to the quantum well planes, 
could be enhanced greatly in a quantum well superlattice 
geometry as a result of an effectively two-dimensional den- 
sity of states for the electrons. They calculated ZT for model 
Bi,Te, superlattices with varying well widths in which the 
potential barriers of the wells were taken to be infinite and 
the widths of the barriers were neglected, and they found 
that ZT increased monotonically with decreasing well width 
attaining values much greater than that for bulk Bi,Te,. They 
also argued that quantum wire superlattices would give fur- 
ther enhancements in the figure of merit [13]. 

The proposals of Hicks and Dresselhaus attracted con- 
siderable attention because of the possibility of obtaining 
greatly improved materials for thermoelectric applications 
from superlattices. However, because these authors did not 
consider the effects of the barrier layers, their calculations 
neglected two physical effects, which are essential to describe 
the figure of merit of realistic superlattice systems even at a 
qualitative level. The first of these effects is that thermal cur- 
rent flows through the barrier layers in addition to the well 
layers, thus producing a decreased ZT. Secondly, electron 
tunneling occurs through the barriers, and the resulting 
superlattice dispersion leads to a decreased density of states 
and correspondingly reduced ZT. 

In order to evaluate realistic composite systems for 
thermoelectric applications, it is necessary to have a descrip- 
tion of their thermoelectric transport properties. This should 
include realistic results for the dependencies of the transport 
properties on the parameters of the systems. Calculations of 
the thermopower [14] and of other properties have been re- 
ported over the years. Recently, several workers have pro- 
vided a realistic description of the figure of merit for model 
superlattice systems [15-19]. 


Quantitative Calculations for 
Superlattice Thermoelectrics 


Lin-Chung and Reinecke [15] included the effects of 
thermal currents in the barrier materials for quantum well 
superlattice systems, which have quasi-two dimensional elec- 
tronic properties in the quantum well regions. They consid- 


10 Naval Research Reviews 


ered Bi,Te,/PbSnTe superlattices and derived the effective 
transport properties as functions of barrier and well widths. 
They found that the thermal currents in the barriers reduced 
the values of ZT at a given well width but that the composite 
superlattice systems nevertheless could have values of ZT 
considerably greater than for the bulk. In this work they did 
not allow carrier tunneling into the barrier regions. Mahan 
and Lyon [20] have also discussed the role of heat transfer in 
the barriers of superlattices. 

Broido and Reinecke [16] included both the effect of 
carrier tunneling for finite potential barriers as well as the 
effect of thermal current flows in the barriers having finite 
widths in order to address these two issues in a unified way. 
A Bi,Te, superlattice was considered in which a fixed poten- 
tial offset V, was present at the boundaries between the wells 
and barriers. The electronic band structure was given by a 
Kronig-Penney description within the effective mass approxi- 
mation, and the transport properties were treated within a 
constant relaxation time approximation to the Boltzmann 
equation using bulk relaxation times. 

They took the superlattice axis to be in the direction of 
the c-axis of the hexagonal unit cell and the x-direction to be 
in the quantum well plane along the a-axis. For current trans- 
port along the quantum well perpendicular to the superlattice 
direction, Fig. 5 illustrates the calculated figure of merit of 
the superlattice, Z.T, as a function of the superlattice period 
d for several different potential barrier heights. Fig. 6 shows 
Z,T as a function of the quantum well width for several ra- 
tios of the barrier width to the well width. In these figures, 
ZT is shown scaled to the corresponding bulk value, which 
is calculated to be Z,, T=0.53. 

From these figures, it is seen that ZT has a maximum 





Figure 5. 


The figure of merit ZT of Bi,Te, superlattice (scaled by Z,,T) 
as a function of the superlattice period d for several potential 
barrier heights (b=a). 























for decreasing well widths. This arises from the effects of 
carrier tunneling, which become greater at smaller well 
widths. For sufficiently small well widths the structure goes 
over into an effectively three dimensional system. From Fig. 
5 one sees that the enhancement in ZT occurs as a result of 
the quantum confinement of carriers, which is greater for 
larger values of V.. For larger V,,, the carrier dispersion along 
the superlattice direction is flatter and the effective density 
of states is enhanced by the confinement. It is seen that Z.T 
can be enhanced by some 50% over Z,,T for modest barrier 
heights of 1 eV. Barrier heights on this order can be obtained 
in semiconductor systems such as InAs/AISb, as well as other 
composite systems involving wider band gap materials such 
as in systems having insulators for the barriers. From Fig. 6 
it is seen that for a fixed V,, the well width at which the maxi- 
mum of ZT occurs depends on the barrier width. This well 
width dependence results from the combined effects of ther- 
mal currents through the barriers and carrier tunneling 
through the barriers. Recently, Sofo and Mahan [19] have 
made a similar treatment of thermoelectric transport in 
superlattice systems. 

The effects of superlattice geometry on semi-metals 
having two carrier components (electrons and holes) such as 
Bi have also been considered [17,21]. A homogeneous sys- 
tem having one type of carrier component has a higher ZT 
than a corresponding two component system because the 
contributions of carriers with different charges tend to can- 
cel in the thermoelectric power. The effects of the superlattice 
potential and of quantum confinement split the electron and 
hole energies into sub-bands; this can transform the material 
into an effective one-carrier system. For example, in 
superlattices such as the Bi/PbTe, the Bi may have a higher 





Figure 6. 


The figure of merit ZT of Bi,Te, superlattice (scaled by Z,,,T) 
as a function of the quantum well width a for several barrier 
widths b for a potential off-set of 200 meV. 




















ZT than in its bulk, semi-metallic form. Such superlattice 
systems involving Bi have been studied [17] using the tech- 
niques developed earlier [16]. A substantial enhancement of 
ZT in the one component superlattice geometry over that of 
bulk Bi was found. It is particularly noteworthy that in this 
case the maximum of ZT occurs for relatively wide wells of 
the order of 60 - 80 A. 


Exploration of Further 
Geometries 


The thermoelectric properties of quantum wire 
superlattices have also been studied [18]. Quantum wire struc- 
tures are now being fabricated in a number of ways, such as 
lithographic techniques, in situ growth of materials in glass 
arrays, and injection of materials into chemical templates. 
Quantum wire superlattices have been studied [18] using tech- 
niques similar to those used to study quantum well 
superlattices [16]. Results for ZT for a rectangular array of 
quantum wires of Bi,Te, for transport along the wire direc- 
tion as functions of potential off-set and quantum wire width 
are similar to those shown for quantum well superlattices 
shown in Figs. 5 and 6. It has been found that the degree of 
enhancement of ZT over that for the bulk is greater in the 
quantum wire superlattices than in the quantum well 
superlattices for a given potential offset above 1 eV. This is 
attributed to the effects of greater confinement in the quan- 
tum wire case [18]. 

Thermoelectric transport in superlattices along the 
superlattice direction has also been considered [22,23]. Re- 
cent preliminary measurements for GaAs/AIAs superlattices 
suggest that the thermal conductivity for transport in this 
direction can be reduced by as much as a factor of ten [24]. 
The physical origin of this reduction is not yet known in 
detail, but it most likely involves scattering due to disorder 
or interfaces. Calculations of ZT of superlattices involving 
Bi,Te, for transport in this direction have been made for a 
Bi,Te, superlattice with several values of the thermal con- 
ductivity, and it has been found that ZT is increased by more 
than a factor of two for decreases in the thermal conductiv- 
ity of a factor of ten. 


Tuning Properties for 
Performance in Composite 
Systems 


Superlattice and other composite systems are attrac- 
tive for thermoelectric applications because a wide variety 
of parameters can be varied to achieve desired properties. 
From the theoretical work to date on quantum well and quan- 
tum wire superlattices we are developing a basic knowledge 
of the thermoelectric transport properties of realistic systems. 
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In particular, it has been found that it is essential to include 
carrier tunneling through the barriers and thermal transport 
in the barriers to describe realistic systems. It has been found 
that a considerable enhancement of the figure of merit, ZT, 
over the corresponding bulk values can be achieved through 
carrier confinement in these systems. This theoretical work 
should stimulate advances in fabricating systems of interest 
as thermoelectric materials as well as further experimental 
work on them. 

There are several important physical effects that need 
to be addressed in describing the thermoelectric transport 
properties of confined systems. It is known that the carrier 
relaxation rates in confined structures depend on the struc- 
ture involved and have complex energy dependencies [25]. 
It is important to include these in describing thermoelectric 
transport in realistic systems. Another important aspect of 
the physics both of confined systems, and also of bulk sys- 
tems, which is not yet fully understood is their thermal trans- 
port. The effects of disorder and interface scattering in these 
systems are expected to decrease thermal transport by 
phonons, and that should produce an improved thermoelec- 
tric figure of merit. Measurements of the thermal conductiv- 
ity in AlAs/GaAs superlattices [26] indicate that it is less 
than either of the bulk values. A theoretical framework for 
understanding these effects on thermoelectric transport is 
needed. In particular, it will be necessary to know how the 
thermal transport in these systems depends on their materi- 
als and structure in order to predict how best to exploit them 
in producing better thermoelectric materials. In summary, 
the basis of a realistic theoretical treatment of thermoelec- 
tric transport in superlattice systems has been developed, but 
key physical issues still need to be investigated. 


Outlook 


Exploration of the huge parameter spaces required to 
find specific materials and structures that provide sought af- 
ter properties remains one of the most daunting challenges 
in materials research. Nonetheless, there has been a great 
deal of progress, both on the experimental front, where new 
or greatly improved tools for controlled synthesis and char- 
acterization have been developed (e.g. molecular beam epi- 
taxy for growing superlattices), and, as discussed in this ar- 
ticle, on the theoretical and computational side. These ad- 
vances are beginning to allow a new paradigm in which ad- 
vances in materials research come from a strong synergism 
between theory, computation, characterization and synthe- 
sis efforts. As discussed in this article, theory has been used 
to advantage in understanding properties of skutterudite based 
materials and in understanding the TE properties of 
superlattice and composite systems. Ongoing efforts include 
the application of first principles methods to other novel 
materials with TE potential, such as the newly synthesized 
Cu-Te compounds discussed by Kanatzidis in this volume 
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and detailed predictions of TE properties in novel 
nanostructured geometries in conjunction with experiment. 
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Thermoelectric 


Materials: 


Solid State 


Synthesis 


Mercouri G. Kanatzidis, Department of Chemistry, Michigan State University and 
Francis J. DiSalvo, Department of Chemistry, Cornell University 


Introduction 


It is remarkable that during the last 35 years the effi- 
ciency of the best thermoelectric materials has not improved 
to any significant extent. This is despite spectacular advances 
in some areas of materials research and technology in that 
period: semiconductor devices, advanced ceramics, super- 
conductors, superalloys, laser, optical and non-linear optical 
materials, photovoltaics and composites, to name but a few. 
The efficiency of thermoelectric coolers operating near room 
temperature is only about 10% of Carnot efficiency; yet the 
thermodynamics of thermoelectric cooling suggest that 
achieving close to 100% of Carnot efficiency may be pos- 
sible. To put that in context, modern freon-based home re- 
frigerators operate at 30 to 35% of Carnot efficiency and 
large building air conditioning units at over 80%. 

An efficient thermoelectric device is fabricated from 
two materials, one n-type and the other a p-type conductor. 
Each material is separately chosen to optimize the figure of 
merit, z, where z = S*o/x; S is the thermopower, 6 the elec- 
trical conductivity, and « the thermal conductivity (1). All 
three of these materials properties are determined by the de- 
tails of the electronic structure and scattering of charge car- 
riers (electrons or holes) and thus are not independently con- 
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trollable parameters. K also has a contribution from lattice 
vibrations, K,,, the phonon thermal conductivity. Thus K = 
K, + K,,, where kK, is the carrier thermal conductivity. 

In order to carry a heat flux of reasonable magnitude, 
moderate to high carrier densities are needed (small band 
gap semiconductors with carrier densities of 10'*/cm?’ to 
metals at 10°?/cm’). In this range of carrier densities, K,, usu- 
ally is much larger than « oh” In that case, the ratio K/o is 
given approximately by the Wiedemann-Franz law, K/o = 
L.T, where L, has the value 2.45 x 10 * watt-Ohm/K? (2). 
Indeed in the best thermoelectric materials currently known, 
Bi,Te, and PbTe, the Wiedemann-Franz law is almost ex- 
actly obeyed in the range of temperatures of optimal device 
performance. For the Bi,Te, material used in commercial 
devices (3), the value of zT at room temperature is about 0.9 
(T is the absolute temperature; z always appears in combi- 
nation with T in engineering calculations.). As zT increases, 
the cooling efficiency increases, a little less than linearly in 
this range of zT. 

Thus improving device performance means improv- 
ing ZT = S*/L,,, or increasing S while keeping moderate to 
large carrier densities of one carrier type. Since the 





thermopower of optimally n or p doped Bi,Te, is on the or- 
der of + 220 'V/K (3), significant improvements in thermo- 
electric cooling efficiencies - to 20 or 30% - could occur if 
materials of reasonable carrier density with thermopowers 
of 300 to 350 p.V/K can be found. 

In this article we focus on the efforts to synthesize bulk 
materials with figures of merit two to four times higher than 
those attainable with Bi,Te,. There are two basic approaches: 
the synthesis of new, structurally and/or compositionally 
complex, materials and the exploitation of new phenomena 
in solids that produce high thermopowers. Both approaches 
rely on advanced synthetic capabilities, on the now routine 
capability to determine complex crystal structures (very dif- 
ficult and extremely time consuming 35 years ago - comput- 
ers have made it all possible), and on automated measure- 
ment of electrical and thermal properties. 

Along with chemical and structural complexity usu- 
ally comes the degrees of freedom to produce and control 
more complex electronic structures. This is certainly evident 
in the field of high temperature superconductors, where the 
structures and six or more element compositions of the ma- 
terials are well adapted to chemical control and modifica- 
tion of the properties. Boltzmann transport theory describes 
both electrical and thermal transport in the vast majority of 
solids (4). This theory provides a general understanding of 
the thermopower that is expressed in the following equa- 
tion: 


S = n’k,’T/3e[d(Ino(E)/dE] at E = E, 


O(E) is the electrical conductivity determined as a function 
of band filling or Fermi energy, E,.. If the electronic scatter- 
ing is independent of energy, then O(E) is just proportional 
to the density of states at E. In the general case S is a mea- 
sure of the difference in 6(E) above and below the the Fermi 
surface - specifically through the logarithmic derivative of o 
with E. 

Since the thermopower of a material is a measure of 
the asymmetry in electronic structure and scattering rates 
near the Fermi level, we desire to produce complexities in 
either or both in a small energy interval (a few kT) near E,. 
While simple materials usually have simple band structures, 
in a few materials interactions between conduction electrons 
or interactions of conduction electrons with local magnetic 
spins lead to sharp features in the electronic excitations near 
the Fermi level that sometimes produce high thermopowers. 

Here, then, we discuss these two approaches to new 
thermoelectric materials, exploring complex materials and 
investigating relevant electron interaction phenomena, by 
choosing examples from our current research. 


Synthetic Methodologies 


a. Complex Materials from the 
alkali-metal polychalcogenide 
flux method 


Structural and compositional complexity can result in 
corresponding complexities in the electronic structure. Such 
complexities may produce the required asymmetry in o(E) 
to obtain large thermopowers. The phonon contribution to 
the thermal conductivity can also be lowered by such struc- 
tural complexity as well as by choosing heavy elements as 
components of the material and by choosing combinations 
of elements that normally make moderate to weak chemical 
bonds. These latter two factors lower the Debye temperature 
and thereby decrease the thermal conductivity. This is just 
what happens in Bi,Te, and PbTe. In general, these latter 
conditions will be met by compounds of main group ele- 
ments beyond the third period. Since the known materials 
that are used in thermoelectric devices are chalcogenide com- 
pounds (the chalcogenides are sulfur, selenium and tellu- 
rium), we are searching for more complex semiconducting 
compounds of this type. 

Since many chalcogenides - especially those that are 
complex - decompose at relatively low temperature (perhaps 
500 to 1000 °C), a low temperature synthesis method is 
needed. We have been exploiting polychalcogenide fluxes 
for exactly these purposes. 

The use of molten salts as solvents for reaction and 
crystal growth of solid state materials has expanded in the 
last few decades to a highly versatile synthetic method. Such 
media have been employed for well over 100 years for high 
temperature single crystal growth but only recently for new 
materials synthesis (5). Although many salts melt at high 
temperatures (2 700 °C), eutectic combinations of binary salts 
and salts of polyatomic species often have melting points 
well below the temperatures of classical solid state synthe- 
ses, making possible their use in the exploration of new chem- 
istry at intermediate temperatures. In some cases, such salts 
act not only as solvents, but also as reactants, providing spe- 
cies which can be incorporated into the final product. 


a.1. Flux Choice 


One class of low melting salts that are especially in- 
triguing is the alkali-metal polychalcogenides, A,Q. 
(A=alkali metal, Q=chalcogenide). Both the solids of these 
compositions and the melts contain Q ” anions, in which x 
Q atoms form a covalently bonded one dimensional chain. 
As solvents for intermediate temperature reactions, A,Q_ salts 
are especially well suited because the melting points range 
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between the extremes of K,S at about 850 °C to K|S, at 145 
°C, with the majority of compositions melting at <300 °C. 
Although few detailed phase diagram studies of the other 
congeners have been performed, the melting points of many 
polychalcogenide salts are lower than 550 ° C. The mini- 
mum melting points are higher for the heavier chalcogenides, 
being about 250 °C for the polyselenides and 350 °C for the 
polytellurides. Low melting, A,Q, fluxes remain nonvola- 
tile over a wide temperature range, and so once above the 
melting point, reaction temperatures can be varied consider- 
ably without concern for solvent loss. 

Reactions between metals and molten A,Q_ are per- 
formed in situ. The powdered reagents (polychalcogenide 
and metal or metal chalcogenide) are mixed under inert at- 
mosphere and loaded into reaction vessels of either pyrex or 
quartz. Once evacuated, the tubes are flame sealed and sub- 
jected to the desired heating program in a computer con- 
trolled furnace. 

The polychalcogenide salts can be made by reacting 
stoichiometric amounts of elemental chalcogenide with al- 
kali metal in liquid NH, at -70 °C. To synthesize new com- 
pounds, one or more metals are added directly to the A,.Q/ 
Q reaction mixture and heated in a sealed pyrex or quartz 
container. Crystalline products either precipitate from the melt 
or form on slow cooling of the melt, depending on the spe- 
cific stoichiometric and processing conditions. 

Polychalcogenide fluxes are highly reactive towards 
metals because they are very strong oxidants. This phenom- 
enon arises from the fact that, although the terminal atoms 
of the (Q.)* chains have a formal oxidation state of 1-, the 
internal atoms are all zero valent and thus are good oxidiz- 
ing agents. Metals are drawn into solution by being oxi- 
dized by the internal atoms of the (Qy chains, which, upon 
being reduced, split into smaller fragments. The solvated 
metal cations are then coordinated by the basic terminal at- 
oms of the (Q )* ligands, forming some manner of soluble 
intermediate, which is followed by nucleation into solid crys- 
tals. Presumably, the nucleated species are in equilibrium 
with the soluble intermediates, especially if the flux is present 
in excess, and hence a solvation/reprecipitation effect (often 
referred to as the mineralizer effect) occurs. This aids in the 
growth of high quality single crystals because the flux can 
redissolve small or poorly formed crystallites and then 
reprecipitate the species onto larger, well-formed crystals. 

At this stage of development of this method, each ter- 
nary metal explored must be treated as a new and unique 
system, and similarities to the reactivity of other metals can 
only be drawn afterwards. Of course, the complexity of the 
resulting compounds may be further increased by adding 
several metals at once to the melt to produce quaternary, 
quintenary, etc. phases. 


a.2. Novel Bismuth Chalcogenides 


Despite the empirical guidelines we have, it is still 
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challenging to choose the particular system for exploration. 
The fact that Bi,Te, is the best material known to date sug- 
gests that it combines many of the necessary features for 
high figure of merit. If there is something special about bis- 
muth in giving rise to simultaneously high electrical con- 
ductivity and thermoelectric power, it should be manifested 
in other compounds of bismuth as well. This suggests a re- 
search direction where we could explore other, more com- 
plex, chalcogenides of bismuth in the hope that some (or all) 
of the key properties would be superior to those of Bi,Te,. 
Further, structurally and compositionally more complex bis- 
muth chalcogenides would, most likely, have a low lattice 
thermal conductivity. This is because a structure with a large 
unit cell is expected for complex materials, which would 
decrease the acoustic phonon mode velocities that are re- 
sponsible for the transfer of heat in materials. The relatively 
weak Bi-Te bonding and the large atomic masses contribute 
as well to the low phonon velocities. Therefore, exploratory 
synthesis in this region of the periodic table becomes a rea- 
sonable activity and, as the preliminary results show, quite 
promising. 

We have obtained many interesting new bismuth con- 
taining phases by reacting the metal with alkali metal chal- 
cogenide salts in polychalcogenide melts. One such phase 
is the orthorhombic BaBiTe,, which was prepared at 600 °C 
from a K,Te, flux in which a mixture of Bi and BaTe were 
dissolved (6). The structure of this material is two-dimen- 
sional with [BiTe,}* layers alternating with Ba cations (Fig- 
ure 1). Many interesting features exist in the structure, one 
of which is the presence of Te-Te bonds between Te atoms 
arranged in ribbons. The flat Te ribbons contain zig-zag Te 
chains and they are alternating between Bi-Te blocks, which 
look like they have been excised out of a rock-salt structure. 





Figure 1. 


The layered structure of BaBiTe,. The two distinct sections 
of the lattice are shown. 














Figure 2. 


Structure of KBi, ,,S,,.. The parts of the structure related to 
Bi,Te, and Cdl, structures are shown. 
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BaBiTe, is a semiconductor with a narrow band-gap 
of 0.35 eV and quite promising electrical properties. The 
electrical conductivity is reasonably high at ~100 S/cm and 
the thermoelectric power reaches +220 'V/K at room tem- 
perature. Therefore, the compound is p-doped. The thermal 
conductivity, which is crucial in assessing the compound’s 
potential, is only 65-70% that of the rhombohedral Bi,Te,. 
The thermal conductivity is suppressed because of the low 
symmetry and because unit cell is much larger than that of 
Bi,Te,. The compound also contains additional soft phonon 
modes from the Ba—Te interactions. The latter are absent 
in Bi,Te,. This result suggests that generally ternary or qua- 
ternary compounds may possess lower thermal conductivity 
than the corresponding binary compounds. In order to fully 
assess the potential of this material for cooling devices, con- 
trolled n and p doping to moderately high levels must be 
achieved. 

Two other interesting phases discovered recently are 
the sulfides KBi,,,S,, and K,Bi,S,, (7). They belong to the 
general family of compounds (A,Q) (Bi,Q,),, (A=alkali 
metal; Q=S, Se) with n=1 and m=6.33, 4, respectively. In an 
architectural context, both compounds can be thought of as 
an intimate composite of two different structure types inter- 
connected to form a 3-D network. They have three-dimen- 
sional orthorhombic structures made up of Bi,Te,-type (NaCl- 
type) blocks and CdI,-type fragments that connect to form 
tunnels filled with eight-coordinate K* cations. This may be 
beneficial to the electronic properties of the compounds which 
may bear similarities to those of Bi,Te,. The lower symme- 
try they possess may result in low thermal conductivity and 
could give rise to a superior thermoelectric material. The 
[Bi, ,,S,,] framework is made of edge-sharing BiS, octahe- 
dra, as shown in Figure 2. The structure of K,Bi,S, , is simi- 





lar but the Bi,Te,-type (NaCl-type) blocks and Cdl. -type frag- 
ments are arranged differently. 

These ternary bismuth sulfides have promising elec- 
trical properties with maximum conductivity and 
thermopower of 200 S/cm and ~90p1V/K, respectively. These 
are unoptimized values, and we believe they can be greatly 
improved by further processing. Using the measured values 
of the electrical resistivity in conjunction with the 
Wiedemann-Franz law, we estimate the maximum possible 
value of the electronic thermal conductivity contribution to 
be below 1% of the total thermal conductivity. Thus, essen- 
tially all heat in these compounds is carried by lattice 
phonons. The pleasant surprise in these compounds comes 
from the rather low thermal conductivity they possess, as 
found especially in the sulfide KBi, ,,S_. Taking as a bench 
mark the room temperature value of the total thermal con- 
ductivity of Bi,Te, (Kk = 1.7 W/m-K), we note that the total 
thermal conductivity of KBi, ,,S,, (8) is actually lower (Fig- 
ure 3). Hence, at least from the perspective of thermal trans- 
port, these compounds satisfy one of the key requirements 
for a useful thermoelectric material, they possess low lattice 
thermal conductivity. This is an important finding because 
KBi, ,,S,,, being a sulfide, is expected to possess higher ther- 
mal conductivity compared to the heavier tellurides. If con- 
trolled doping can enhance the electrical conductivity and at 
the same time preserve or even increase the thermopower, in 
the case of KBi, ,,S,,, we indeed might have a promising ther- 
moelectric material. To achieve this, we need additional in- 
formation regarding the transport properties including car- 
rier concentrations and mobilities. This work is currently 
underway (8). Both KBi,,,S,, and K,Bi,S,, melt with no 





Figure 3. 


The thermal conductivity of polycrystalline cylindrical samples 
of KBi, ,,S,, (open squares)and K,Bi,S,, (open circles) as a 
function of temperature. The arrow indicates the room 


temperature value of the lattice thermal conductivity of Bi, Te.. 
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decomposition at 710 °C and 713 °C, suggesting they will 
be amenable to thermoelectric element fabrication and pro- 
cessing similar to that currently used in Bi,Te, technology. 

While we expect to find many new bismuth chalco- 
genide compounds by this method, it is not necessary to limit 
the search to bismuth. Other heavy main group metals may 
also lead to high thermopower materials. The synthetic search 
is well underway; although much remains to be done. 


b. Electron Interactions: 
Intermediate Valence in Rare 
Earth Intermetallic Compounds 


Many new phenomena have been discovered in sol- 
ids since the early 60’s. Of particular interest is the unusual 
behavior of rare earth atoms in some intermetallic com- 
pounds, since this behavior can include very high 
thermopowers relative to other metallic systems. 

Normally, the f electrons of rare earth atoms are well 
shielded from bonding interactions by the s and d valence 
electrons. The f state of rare earth atoms in compounds is 
only weakly perturbed from that found in free rare earth at- 
oms or ions. In some materials, the f electrons of the rare 
earth elements Ce, Sm, Eu, Tm and Yb interact with the con- 
duction electrons to produce a variety of electronic states 
(9): Heavy Fermions, Kondo Insulators, Kondo Lattices, and 
Intermediate Valence (IV). In the IV state the occupancy of 
the f state appears to be non-integral. While the details of the 
unusual states that result from this conduction electron - f 
State interaction are still being examined by experimental- 
ists and theorists, it is experimentally apparent that these IV 
compounds have very unusual electronic transport proper- 
ties, including large thermopowers. 

CePd, is a metallic conductor and is the prototypical 
IV material with a large thermopower. YbAI, (10) another 
IV material has a thermopower of opposite sign to that of 
CePd, (YbAI, is n-type while CePd, is p-type). Thus IV 
materials exist with both types of carriers, a necessary con- 
dition to obtain devices with the highest cooling efficien- 
cies. The thermopowers of these IV materials are higher than 
for any other kind of metallic system. Note also that the 
thermopower of both CePd, and YbAI, is rather insensitive 
to temperature above 100 K, in contrast to that of metals or 
of Bi,Te, (11) where SaT. Since these [TV compounds are 
metallic, we expect that the ratio K/o will be near the 
Wiedemann-Franz limit. In CePd, this ratio is quite sample 
dependent, probably due to slight differences in stoichiom- 
etry or to the presence of anti-site and other defects. Never- 
theless, the ratio we measure is about a factor of two larger 
than expected. 

The IV state produces considerable changes in the elec- 
tronic and magnetic properties of the material. In most inter- 
metallic compounds the valence electron configuration of 
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Ce is the same as in a free Ce atom: 4f'5d'6s*. The one f 
electron produces a magnetic moment of 2.54 1, (Bohr mag- 
netons) per Ce atom, while the remaining s and d states be- 
come delocalized as part of the electronic band structure. 
Consequently, most intermetallic Ce compounds are para- 
magnetic with a moment that is also typical of Ce*? ions in 
an ionic solid (a 4f' electron configuration for the single re- 
maining valence electron). Thus, even in intermetallic com- 
pounds, when the normal 2.54 1, moment is observed, the 
Ce atom is referred to as trivalent. In a few Ce compounds 
the f moment disappears, probably because of a change in 
the Ce valence electron configuration to 5d*6s?. While such 
a change in configuration costs energy, increased overlap of 
the d orbitals with conduction band states may result, lower- 
ing the electronic energy sufficiently to make up for the ini- 
tial cost of the configuration change. When Ce has no mo- 
ment in such intermetallic phases, it is called tetravalent, since 
the moment is also lost in ionic Ce** compounds, where no 
valence electrons remain. 

This latter view of tetravalent Ce is not the only way 
that tetravalent-like behavior can be observed. In some Ce 
intermetallics and below a certain temperature, called the 
Kondo temperature, the conduction electrons interact with 
the Ce so as to screen the Ce moment. In compounds where 
no moment is observed, it may be that the Kondo tempera- 
ture is well above room temperature. In typical trivalent Ce 
compounds the Kondo temperature is low or zero. 

The IV state is intermediate between typical trivalent 
and tetravalent states. In the IV state the magnetic moment 
is reduced below 2.54 1, and is often temperature depen- 
dent. The IV compounds have room temperature suscepti- 
bilities only half that expected for Ce with the full f' mo- 
ment, and most of the temperature dependence is lost. At 
low temperatures, current theories would suggest that the 
Ce magnetic moment is effectively cancelled (screened) by 
the conduction electrons through a Kondo resonance. The 
IV state is also signaled by changes in all other properties, 
such as electrical conductivity, thermopower and Ce to near- 
neighbor bond lengths. In some cases, before measuring 
thermopower, it is easier to look at other properties to see if 
the IV state is present, and if not, different properties may 
not need to be measured. 

That is the good news. Now for the challenges. The IV 
state is rather rare. It has been identified in a few dozen in- 
termetallic compounds out of perhaps 3 or 4 thousand com- 
pounds containing one of the five rare earths named above. 
The IV state is also rather sensitive to chemical perturba- 
tions. For example, non-stoichiometry or chemical substitu- 
tion for either Ce or Pd in CePd,, cause large changes in the 
physical properties and supress the IV state (see for example 
refs. 12 and 13). 

Given these facts about the known IV materials, many 
questions remain, such as: Can we find other [TV compounds? 
What is the highest thermopower we can obtain from such a 
material? Can we optimize the figure of merit by chemical 





doping and other crystal-chemical techniques in IV materi- 
als - without destroying the IV state? Since intermetallic 
compounds tend to be brittle, can we find ways to make the 
IV materials less so? 

In order to explore the potential of IV materials to be 
used in higher efficiency thermoelectric devices, it is clear 
from the above questions that the synthesis of new and modi- 
fied materials must play a major role. Here, we focus on our 
synthetic techniques and recent results to illustrate how such 
materials are synthesized, and we show some of their crys- 
tallographic or physical properties. 

Intermetallic compounds are usually prepared from the 
metallic elements by melting at high temperatures (from 1000 
to 2500 °C). We typically use an arc furnace, in which an 
electric arc is used to melt the reactants on a water cooled 
copper hearth and the product is rapidly cooled when the arc 
is turned off. With sufficient care the inclusion of Cu in the 
product can be avoided. In other cases, where more control 
of the reaction temperature is needed, an RF furnace is used. 
An RF furnace couples an RF field generated in a cylindri- 
cal coil to an electrically conducting susceptor (usually made 
from Ta or Mo) inside the coil. The susceptor is also the 
sample container, if the susceptor is non-reactive with the 
sample under the preparation conditions. Otherwise the 
sample is put in some other non-reactive container inside the 
susceptor. Particular cases may require a modified approach. 
For example, Yb has a high vapor pressure at its melting 
point, so its compounds are best synthesized in a sealed con- 
tainer. Mo crucibles that are welded shut and then heated in 
the RF furnace are suitable as both container and susceptor 
in that case. 

Depending upon the reactant composition chosen, the 
molten product will freeze either into a single intermetallic 
compound (a congruently melting solid) or the melt will 
freeze into a mixture of solid phases. In most cases anneal- 
ing at modest temperatures (700 to 1000 °C) is needed to 
produce a well-ordered crystalline product and/or the equi- 
librium phases that result from that starting composition. 
When searching for new phases, the compound stoichiom- 
etry is generally not known, but often is in the ratio of small 
integers. By making samples of different starting composi- 
tions and comparing powder X-ray diffraction patterns, we 
can identify the compositions at which a single homogeneous 
phase exists. In many cases the structure of the phases pro- 
duced can be identified by powder diffraction. In other cases 
single crystals are needed so that the structure (and thus the 
composition) can be determined by single crystal X-ray dif- 
fraction methods. By heating a homogeneous material at tem- 
peratures just below the point at which some melting occurs, 
or just below where decomposition occurs, small crystals 
often grow within a polycrystalline mass by grain coarsen- 
ing, and the crystals may be separated by lightly crushing 
the sample. 

Our synthetic program currently has two foci: finding 
new ternary or quaternary rare earth compounds that are [IV 


and modifying known IV materials. We next discuss examples 
of each to give a a better sense of the synthetic challenges 
and the spectrum of materials that we are synthesizing. 


b.1. New Ce Intermetallic 
Compounds 


CePd, and CeSn, are both IV materials. What happens 
if all three elements are included in a ternary compound? 
Will more IV materials be produced? Will mixing the non- 
rare earth elements, especially in an ordered fashion, pro- 
duce larger thermopowers? These kinds of questions lead us 
to the study of the Ce-Pd-Sn phase diagram. The phase dia- 
gram, when completely explored, indicates the stoichiom- 
etry and structure of the compounds that exist and often in- 
dicates other information such as melting points, etc. All the 
binary compounds of Ce-Pd and Ce-Sn were already known, 
as well as two ternary phases: CePdSn and CePd,Sn,,. In our 
studies we discovered several new ternaries, but discuss only 
one here. 

When a small amount of Sn is added to CePd,, an un- 
usual compound results, Ce,Pd, Sn (14), which crystallizes 
in a new structure type. The structure is an ordered and partly 
filled superstructure of CePd,. This compound is surprising 
in that it is uncommon to find intermetallic phases with sto- 
ichiometric coefficients so much larger than 1. The compound 
is produced by arc melting and annealing at 900 °C for sev- 
eral weeks. Interestingly, we have discovered that Sn may 
be replaced by a wide variety of main group metals. So far, 
the resistivity of two compounds, Ce,Pd,,Ga and Ce,Pd, In, 
indicate strong interactions with the conduction electrons 
(figure 4). Magnetic studies and thermopower measurements 





Figure 4. 


The resistivity of Ce,Pd,,In and Ce,Pd,,Ga show considerable 
increases in resistivity over that of La,Pd,,In, a compound 
with identical structure but no f electrons. It appears that the 
Ce compounds have Kondo temperatures on the order of 20 
K, but since the room temperature thermopowers are low 
and Ce magnetic moments are normal, they do not qualify 
as IV materials. 
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show, however, that the strong IV character of CePd, is lost 
with even this small perturbation on the composition and 
structure. 

Even though we have synthesized several dozen new 
Ce intermetallics, we have so far found hints of only one 
new IV compound, Ce,CoSi, (15). At present there are no 
theoretical or even empirical guidelines which help in de- 
ciding which particular materials will be IV, and, given the 
small number of known IV materials in the family of Ce 
compounds, the expected success rate of discovering new 
ones is low. The magnetic susceptibility of the series 
Ce,Au, Co Si, (figure 5) shows the evolution of the behav- 
ior from typical Ce** moments at small x to apparent reduced 
moments at x = 1.0. Ce,CoSi, is quite brittle in polycrystal- 
line form, and we have not yet been able to determine the 
thermopower. Since single crystals should be mechanically 
more robust, we are attempting to prepare them for 
thermopower and other property measurements. Also since 
Ce,CoSi, has a hexagonal structure, its properties are ex- 
pected to be anisotropic, and single crysta! measurements 
will be necessary to investigate such anisotropy as well. 


b.2. Modifying known IV materials 


We are also trying to understand the doping behavior 


of the known IV materials to see how the thermopower is 
affected by such chemical manipulation. In particular, some 
theoretical models suggest that the thermopower should be 
independent of Ce concentration in an IV compound, if all 

j i . Our 
studies of the dilution of Ce in CePd, show that the situation 
is complex and not yet completely understood. Figure 6 
shows the room temperature thermopower of several differ- 
ent alloy systems: Ce, Y Pd, (16), Ce, (La, ste. (17), 
and Ce, Nd Pd, (17). To prepare homogeneous samples of 
such complex alloys, great care must be taken. In many cases 
reported in the literature, such care is clearly not taken and 
essentially meaningless studies are performed on 
inhomogeneous materials. When such alloys are melted they 
always form with composition gradients on solidifying. These 
gradients can be rather large. The composition gradients can 
only be removed after melting by prolonged annealing at 
temperatures approaching the point at which some melting 
first occurs. At this point the diffusion coefficient is usually 
large enough that, for gram size samples, the material be- 
comes homogeneous in a few days. Such synthesis strate- 
gies will not work for large samples. On a technological scale 
a powder metallurgical approach could work, if processing 
temperatures are low enough to prevent the formation of any 
melt. 





Figure 5. 


The magnetic susceptibility of Ce Au, Co Si, shows typical IV like susceptibility at x = 1.0. At the Au rich end of the alloys 


(small x), typical trivalent behavior is observed. 
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Figure 6. 

The thermopower of Ce, Y Pd, (@), Ce, (La, .Y,.),Pd, (+), 
and Ce, Nd,Pd, (4) is shown vs. x. The solid line is the 
theoretcial dependence based on the single Ce model and 
the Nordheim-Gorter rule. 
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The data forCe, Y Pd,, taken from the literature, show 
a decline of the thermopower with doping. Such doping also 
leads to a decrease in the unit cell volume. The decrease in 
thermopower has been interpreted as due to the chemical 
“pressure” caused by the replacement of Ce by the smaller 
Y. Pressure favors the formation of tetravalent Ce, since it is 
also smaller than trivalent Ce. 

If chemical pressure were the only reason for the re- 
duced thermopower, then replacement of Ce by the correct 
mixture of La (larger than Ce) and Y (smaller than Ce) would 
keep the crystallographic unit cell exactly the same size. The 
data on such an alloy shows that the thermopower decreases 
nonetheless, but not quite as rapidly as for Y alone. This re- 
duction, however, might be due to atomic scale fluctuations 
in the lattice strain due to the different size of La and Y. Such 
fluctuations would not be apparent in the normal x-ray dif- 
fraction pattern where the lattice parameters are found to be 
independent of x. 

In order to avoid lattice strain effects, we looked at 
dilution of Ce with Nd. NdPd, and CePd, have lattice pa- 
rameters that differ by a very small amount (the cubic lattice 
parameters for both compounds are 4.128(2) A). However, 
the thermopower again decreases with increasing Nd con- 
centration, but not quite as rapidly as for the previous two 
materials. It would appear then that the Ce electronic state is 
not independent of concentration when we control the lat- 
tice and electronic features to be as close to those of CePd, 
as is possible. This may indicate that the effects we seek to 
control are collective, and not single Ce ion effects. How- 


ever, further studies and comparisons with theoretical mod- 
els are needed in order to fully understand the implications 
of these measurements. 

Our exploration of IV materials as potential thermo- 
electric devices is just beginning. Our fundamental under- 
standing of the IV phenomenon is elementary; our ability to 
predict the composition of new IV materials is virtually nil; 
and our understanding of how to manipulate the IV state by 
doping is just starting. We have a considerable way to go to 
see if IV phenomena will be useful for the purpose of ther- 
moelectric cooling. 


Summary 


Bulk Bi,Te, is a relatively simple small bandgap 
semiconductor whose electronic properties have been un- 
derstood both experimentally and theoretically for more than 
30 years. It’s thermoelectric performance is near optimal for 
a bulk material with a single conventional conduction and 
valence band. Yet it is clear that we still cannot predict which 
specific materials will have significantly better thermoelec- 
tric properties. However, as discussed in the text, we do have 
general “guidelines” that help focus our attention on broad 
classes of materials. Generally, materials with electronic 
structures that are complex in a narrow range of energy near 
E, are the focus of the searches underway. 

New synthetic methods to prepare complex materials, 
an ability to easily determine the structure of those phases, 
and recent discoveries of new phenomena in conducting sol- 
ids give us hope that we will someday replace the Bi,Te, 
family as the highest efficiency thermoelectric materials. We 
as a materials community have yet to find the desired mate- 
rials. There are promising leads, but the task is challenging 
and time intensive. The search is so far highly empirical: it is 
guided by our past knowledge, but it is not a predictable or 
easily definable path. We hope that a continued and growing 
interaction with theorists may provide better guidance in this 
process. The potential payoff is enormous in both economic 
and environmental terms; yet the “risk” is high: it may take 
many years of focused research to find such materials. And 
as is usually the case, it will likely take a considerable devel- 
opment and engineering effort before new devices will be 
affordable and generally available. However, without the 
ability to synthesize and characterize high quality materials 
of considerable complexity, our chances of success are very 
slim indeed. 
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The Thermoelectric 
Properties of 
Skutterudites 


Glen A. Slack, Department of Physics, Rensselaer Polytechnic Institute, Troy, New York and 
Jean-Pierre Fleurial and Thierry Caillat, Jet Propulsion Laboratory, Pasadena, California 


Introduction 


Skutterudites are a class of compounds based on CoAs., 
which is actually the mineral skutterudite. This name comes 
from the Norwegian town of Skutterud, near the location 
where the mineral was originally found. There are nine bi- 
nary semiconducting compounds in this group: CoP,, CoAs,, 
CoSb,, RhP,, RhAs,, RhSb,, IrP., IrAs,, and IrSb,. in some 
ways ‘they resemble the III-V semiconductors AIP, AlAs, 
AISb, GaP, GaAs, GaSb, InP, InAs, and InSb. Both sets of 
compounds are covalent semiconductors with cubic crystal 
structures. They can be made into both n-type and p-type 
variants by suitable doping techniques. 

One of the requirements for a good thermoelectric com- 
pound is a high carrier mobility which correlates with a small 
electronegativity difference between the elements in the com- 
pound “”. The smaller the value, the more nearly covalent 
the compound. Skutterudites are significantly more covalent 
than the III-V compounds. A number of recent papers have 
appeared describing the thermoelectric properties of the 
skutterudites °°”). 


Figure of merit 


The performance of a material in a thermoelectric re- 
frigerator or generator can be related to the dimensionless 
figure of merit, ZT, where T is the absolute temperature, and 
Z at the temperature T is computed from: 

Z=S*o/A (1) 
where S= the absolute Seebeck coefficient, o= the electrical 
conductivity, and A= the thermal conductivity. From Eq.(1) 
we can see that a good thermoelectric material will have a 


high electrical conductivity and a low thermal conductivity. 
The Wiedemann-Franz-Lorenz law states that for a metal: 


A, =L,oT (2) 
where A. is the electronic component of the thermal conduc- 
tivity, and L, is the Lorenz number. In a semiconductor the 


lattice phonons or quantized lattice vibrations contribute sig- 
nificantly to the total thermal conductivity so that 
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A(TOTAL) = A, + A, (3) 


Here A= the lattice thermal conductivity. In the case of semi- 
conductors we have: 


71 = S’al_-_S*oT . _S* 
ActAg LooT+Ag Lo+(Ag/oT) = (4) 





Eq.(4) shows the influence of A, on the figure of merit. 

As explained by Slack“ the ideal semiconducting 
material for a thermoelectric device is one which is a “phonon 
glass and electron crystal” or PGEC. A phonon glass has the 
very smallest possible value of A, because the atomic disor- 
der reduces the phonon propagation distances to their abso- 
lute minimum value. Hence A, becomes A,,,.. The electron 
crystal should appear to the electrons (or holes) as a very 
perfectly ordered, covalent solid so that the electrons are not 
scattered very much as they propagate through the crystal. 
This dual nature of disorder for phonons and order for elec- 
trons is the critical feature for finding new and better ther- 
moelectric materials. We shall show that skutterudites have 
a very special nature that may make this possible. 


Low Lattice Thermal 
Conductivities 


The best currently employed thermoelectric materials 
are Ge-Si, Bi,Te,-Sb,Te,, PbTe-PbSe and Bi-Sb mixed crys- 
tals. The mixed crystals have higher ZT values than either 
parent because the atomic mixture produces what is termed 
“mass fluctuation” scattering of the lattice phonons. This 
effect was discovered by Eucken and Kuhn in 1928. This 
randomness of the mass, size and charge of the atoms on a 
particular lattice position in the crystal can also tend to scat- 
ter electrons as well. Thus mixed-crystal thermoelectrics 
have to be carefully tailored to produce some decrease in A 
without too much of a decrease in 6. Mixed crystals can be 
made between the various III-V compounds, and they can 
also be made in the skutterudite system. 

The special structure of the skutterudites offers a sec- 
ond method of reducing A, that is not found in the III-V 
compounds or any of the other currently used compounds. 
Each cubic unit cell of the skutterudite structure contains 
two ordered voids”. These voids can be populated with “ex- 
tra” atoms whose job it is to scatter phonons but not elec- 
trons. The possibility of tailoring these “extra” atoms to yield 
a PGEC behavior of the skutterudites is creating much cur- 
rent interest in these materials. 

Both n-type and p-type conductivity samples have been 
obtained for skutterudites, using several preparation tech- 
niques’. Associated with a low hole effective mass, very 
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high carrier mobilities, low electrical resistivities and mod- 
erate Seebeck coefficients are obtained in p-type 
skutterudites®. For a comparable doping level, the carrier 
mobilities of n-type samples are about an order of magni- 
tude lower than the values achieved on p-type samples. But 
the much larger electron effective masses and Seebeck coef- 
ficients make n-type skutterudites even more promising can- 
didates. However, the thermal conductivities of the binary 
skutterudite compounds are too large (100-150 mW/cm K at 
300K), particularly at low temperatures, to be useful for ther- 
moelectric applications. Substantial reductions in the lattice 
thermal conductivity must be obtained to achieve values com- 
parable to those of state of the art thermoelectric materials 
(10-40 mW/cm K at 300K). Several approaches to a low lat- 
tice thermal conductivity in skutterudites are presented in 
the next section. 


Reducing the Lattice 
Thermal Conductivity 


Because of the importance of A, on the value of ZT a 
great deal of effort is being devoted to understanding the 
mechanisms that may be employed for reducing the A, of 
skutterudites. These may be listed as methods for scattering 
the lattice phonons: 


Mass and Strain Fluctuation 

Void Fillers 

Charge Carriers from Dopants 
Mixed-Valence Cations 

Inner-Shell Excitations (d-shell & f-shell) 
Grain Boundaries, Precipitates, Others 


We shall discuss a number of these possibilities. 

The phonon heat transport in crystals occurs over the 
whole range of phonon energies from zero up to an energy 
of k®, where k = Boltzmann’s constant and 8, = the Debye 
temperature. Each one of the phonon scattering mechanisms 
listed will scatter phonons in some particular energy range 
while interacting weakly with all of the other phonons. For 
example the mass fluctuation scattering cross-section in- 
creases as V where v is the phonon frequency. It is termed 
Rayleigh scattering. Thus it is very effective on phonons with 
energies near k@, but allows the low frequency phonons to 
travel almost unimpeded. We can only reduce A, to Avo if 
phonons of all energies have a mean-free path equal to one- 
half their wavelength. Thus simple mass fluctuation scatter- 
ing, produced by making mixed crystals or alloy crystals, 
can never yield A,,,... The best hope is to combine mass fluc- 
tuation scattering with another mechanism that interacts with 
the lower frequency phonons. This technique was employed 
by Dismukes et al." in Si-Ge samples where the lowest 
frequency phonons were scattered by interaction with the 





free carriers in these degenerately doped materials. The mass 
fluctuation scattering was produced by using samples with 
70 atom % Si and 30 atom % Ge. The combination of these 
two mechanisms produces quite low thermal conductivities, 
but these materials still do not possess A, = A\gy 88 Shown 
by Slack and Hussain™. 


Intrinsic Thermal 
Conductivity 


Before we discuss in detail the results on reducing the 
A, of skutterudites we need to know how the pure crystals 
behave. In Fig. 1 we show the present data on the lattice 
thermal conductivity of CoSb,, RhSb,, and IrSb,. Data on 
the other skutterudites is quite limited. Note in Fig. 1 the 
almost exponential rise in A as the temperature decreases 
below 300K. This feature is characteristic of pure, isotopi- 
cally clean crystals for T<@,. Since 6, is close to 300K for all 
three crystals, this rise starts below room temperature. A simi- 





Figure 1 


The Lattice Thermal Conductivity of Pure, Undoped CoSb, 
and RhSb, crystals and polycrystalline IrSb, Versus 
Temperature. The RhSb, curve stops at 300K, the CoSb, at 
1000K. 
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lar rise is observed in InSb” below its Debye temperature 
of 207K. The decreasing A below 10K in Fig. 1 is caused in 
CoSb, by phonons being scattered from the walls of the 1 
millimeter diameter single crystal. In RhSb, the decrease is 
caused by some unknown internal defects, possibly Sb pre- 
cipitates. In the polycrystalline IrSb, it is caused by phonon 
scattering at the boundaries of the 7 1 diameter grains. In 
Fig. | it is apparent that the A of CoSb, rises above 700K. 
This rise is caused by the bipolar diffusion of carriers which 
are thermally excited across the small band gap (Eg ~ 0.63 
eV) of CoSb,. The band gap of IrSb, is too large for this 
effect to appear. Such bipolar heat transport has been well 
documented” in Si and Ge at high temperatures. 


Minimum Thermal 
Conductivity 


Every crystal has a lower limit to its lattice thermal 
conductivity. This concept and the method of calculating this 
value of A,,,., has been explained by Slack‘. The minimum 
value will be attained when all of the acoustic phonons have 
a mean-free-path equal to one-half of their wavelength and 
all optic phonons have a mean-free-path equal to the aver- 
age interatomic distance. The sum of these two contribu- 
tions gives A,,,.,. The calculated values of A,,,,, for IrSb, ver- 
sus temperature are given in Fig. 2. The values for IrSb, and 
CoSb, are almost identical. 





Figure 2. 


The calculated minimum thermal conductivity versus 
temperature for IrSb.,,. 
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Mass and Strain 
Fluctuation Scattering 
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The mass and strain fluctuation scattering of phonons 
in crystals is usually produced either by isotopic substitu- 
tion or isoelectronic substitution of one or several elements. 
In the skutterudite compounds solid solutions have been re- 
ported to exist in the IrSb, - RhSb, and CoAs, - CoP, sys- 
tems over the whole range of compositions “'*'», Recent 
results at the Jet Propulsion Laboratory show that the IrSb, - 
CoSb, system also forms a complete series of solid solu- 
tions. In the CoSb, - CoAs, system there is a wide miscibil- 
ity gap" over the central 80 mole % of the system. These 
and other recent results''*'” are summarized in Table 1. 

An example of the possible thermal conductivity re- 
duction from mixed crystals is useful. Mixed crystals of 
RhSb, - IrSb, can be made over the whole composition range. 
The thermal conductivity of mixed crystals of 50 mole % 
RhSb, with 50 mole % IrSb, has been measured above room 
temperature. A careful analysis of these results shows that 
the decrease in A is caused only by the mass difference be- 
tween the Rh and Ir, there is no strain-fluctuation compo- 
nent. Thus Rh behaves much like a very light isotope of Ir. A 
similar result has been found in AlAs - GaAs mixed crystals 
by Afromowitz'*’, where Al acts like a light isotope of Ga. 
The Rh, .Ir, Sb, crystal has a thermal conductivity at 300K 
of A,(0-50) = 90 mW/cm K. This is 68% of that of the 
average of the two parents, RhSb, and IrSb,, and is much 
larger than A, of 3.1 mW/cm K. Some recent data on mixed 
crystals at 300K are shown in Figure 3, where measured val- 
ues of the lattice thermal conductivity are given. The mea- 
sured total A has had the electronic component subtracted, 
and also has been corrected for smail porosity effects. Thus 
the numbers correspond to theoretically dense samples. 

The effect of the mass and strain fluctuation scattering 
of phonons by these point defects in skutterudites can be 





Table 1. 


Some Known Skutterudite Mixed Crystals 





Compound | Partial Range of Solid| Full Range of Solid 


CoP, 
CoAs, 
CoSb, 
RhSb, 
IrAs, 


CoAs, 
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calculated using the theory developed by Callaway''®’. This 
theory permits us to calculate A, versus chemical composi- 
tion over the whole composition range. Some of the results 
are shown in Fig. 3 for various CoSb, mixed crystals. Curves 
2 and 3 are for samples that are miscible over the whole 
composition range while the other five curves are predic- 
tions for samples that have a limited range of actual misci- 
bility. We have found, for example, that CoSb, will accept 
up to 12 mole % IrP,, see curve 7. This means that the ex- 
pected A, will be about 13 mW/cm K at 300K. This reduc- 
tion in A, is quite substantial compared to pure CoSb, at 104 
mW/cm K. Further studies of these effects are presently un- 
derway. 





Figure 3. 


The Calculated Lattice Thermal Conductivity of CoSb, Mixed 
Crystals at 300K as a Function of the Mole Percent of CoSb,, 
The circular points have been measured for CoSb.-IrSb,, the 
square point is for the ninth sample (CoSb, -Fe, .Ni, Sb.) in 
Table 7 with 79 mole % CoSb., 
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Scattering by Void 
Fillers 


Numerous different atoms have been introduced into 
the voids in skutterudite crystals”. Among these are Ba, Sr, 
La, Ce, Th, U, Nd and Sm. In IrSb,, for example, Ba** ions 
just fill the voids with no space left over. If all of the voids in 
the crystal, which are always periodically arranged, are filled 
with Ba atoms so that the chemical formula is Ir,BaSb,,, the 
phonon propagation will detect no disorder and the thermal 
conductivity will be high, i.e. about 160 mW/cm K at room 
temperature. If we fill only one-half of the voids with Ba 
and leave the others vacant, than we will have mass fluctua- 
tion scattering similar to that discussed previously. This ef- 





fect scatters predominantly the higher frequency phonons. 
However if we place smaller size atoms in the voids such as 
Nd* or Sm*, these “rattle about” because the voids are sig- 
nificantly larger in diameter than the Sm or Nd ions. Thus 
looser bonding of undersize atoms will change (decrease) 
the frequency of the lattice phonons that are being scattered. 
Furthermore, the undersize atoms will not vibrate in phase 
with one another. This means that all of the voids can be 
filled with undersize atoms and the atomic “disorder” caused 
by the uncoordinated “rattling of the atoms in their cages” 
will scatter the lower frequency lattice phonons. The result 
is a drastic reduction in the thermal conductivity as shown in 
Fig. 4. Here La, Nd, or Sm atoms in the voids in IrSb, have 





Figure 4. 


The Lattice Thermal Conductivity of Pure, Unfilled, 
Polycrystalline IrSb, and of Rare-Earth Filled IrSb, Versus 
Temperature The curve for A,,,,is also shown. See reference 
#20. 
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produced the phonon scattering. The lattice thermal conduc- 
tivity at room temperature has dropped from 160 mW/cm K 
to 11 mW/cm K for the Nd-doped sample, a factor of 15 
decrease. Further research is presently underway to deter- 
mine whether these void occupants can be combined with 
mass fluctuation and strain fluctuation scattering in order to 
further reduce the lattice thermal conductivity. In addition 
other void occupants such as Bi, Ar, Kr, Xe, K, and Ti are 
being considered. 


Charge Carriers from 
Dopants 


The effects of the charge carriers from n-type and p- 
type dopants on the thermal conductivity of semiconductors 
has been studied in Ge, Si, SiC, GaSb, and GaAs by numer- 
ous authors‘!!! In Si, both n and p dopants, when present 
at high concentrations so that the semiconductors are degen- 
erate, are effective as scatterers of the very low frequency 
phonons. In the skutterudites IrSb, and CoSb, we have found 
that n-type dopants in degenerately doped samples are ef- 
fective scatterers of low frequency phonons. Fig. 5 shows 
that at temperatures below 200K the thermal conductivity is 
drastically reduced. P-type dopants seem to produce a much 
smaller effect. The reason for the difference is not currently 
known. 

The charge-carrier scattering of phonons has been 
shown to improve the efficiency of Ge-Si thermoelectric gen- 
erators. Thus it could be useful in the n-type skutterudites 





Figure 5. 


The Lattice Thermal Conductivity of Pure, Unfilled, 
Polycrystalline IrSb, and single crystal CoSb, Doped with 
Various Levels of n-type and p-type Impurities Versus 
Temperature. 
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when coupled with void-filling and mass-fluctuation scat- 
tering. 


Mixed-Valence Cations 


In the simple binary skutterudites like IrSb, the iri- 
dium atoms each donate 3 electrons to the covalent bonding, 
the antimony atoms each donate 5. Thus there are 18 bond- 
ing electrons donated per formula unit. There are nine bonds. 
Each bond is made up of two electrons of opposite spin. This 
same electron counting scheme is applicable to the semi- 
conducting ternary skutterudites. There are a number of ter- 
nary skutterudites that contain Fe, Ru, or Os atoms. These 
atoms tend to have mixed valences similar to iron in Fe,O, 
or ruthenium in RuSb,-RuTe, mixed crystals. In Fe,O, iron 
occurs as both Fe”* and Fe**. In RuSb, the ruthenium i is Ru 
while in the RuTe, it is Ru*, in the mixed crystals both va- 
lences occur. Other examples are the nominal skutterudites 
FeSb,Te, RuSb,Te, OsSb,Te, and Ru, Pd, Sb,, which are 
never exactly stoichiometric and possess metal i ions with two 
different charge states. 

Let us write RuSb,Te more correctly as: 


0) Ru* Ru* [Sb,Te,], 


We find experimentally that some vacancies occur on 
the metal sites, hence the square-denoting a vacancy - witha 
fractional concentration v. The sums v+x+y = | and a+b = 1 
are equal to unity. In one case we have replaced some Sb 
atoms by Ge atoms on the antimony sites so that the compo- 
sition is: 


0) Ru* Ru* [Sb,Te,Ge.]}, 


where: a+b+c = 1. For the skutterudite Ru, <Pd, 5b, we have: 


0) Ru* Ru* Pd*,Sb, 


The compositions are collected in Table 2 along with 
the measured value of the lattice thermal conductivity at room 
temperature. 


Model of Phonon 
Scattering by Mixed- 
Valence Cations 


In this model the phonons are assumed to be scattered 
when an electron pair in the d-shell of Ru** is transferred to 
a neighboring Ru**. So one ruthenium atom goes from 4d° - 
> 4d* while a neighboring one goes from 4d* -> 4d°. In the 
process a phonon is scattered. The scattering depends on the 
presence of both charge states. Thus the phonon scattering 
probability and hence the lattice thermal resistivity, W ’ 
should vary as the product of c(Ru**).c(Ru**) where c is the 
concentration in the atoms per om volume and W_ = I/A,. 
Since each unit cell of volume a,’ contains 8 possible metal 
atom sites, then 


c(Ru**) = sy 
a 


and 


c(Ru2*) = ax 
ap 


and c(Ru**).c(Ru**) = 64xy/ a,° 
tion parameter, P, as 


. Let us define a concentra- 


64xy 
af 


P = 





Table 2. 


Composition Parameters and Thermal Conductivity at 300K for Several Ruthenium Skutterudites. 
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Number 


Nominal A, 
Compo- milliW/ 
sition cm K 





AhSb, | AhSb, 105 





OX27 | RuSb,Te ; 33 





RPS6 |Ru,.Pd,.S 28.5 


b, 





OX51 | AuSb,Te 





OX52 |Ru(Ge-Sb- 
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*P in units of 10°/cm® 
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We have given this in Table 2 in units of (10“/cm*). In 
Fig. 6 we plot W, versus P. At P=0 we have employed the 
measured value of W, for RhSb, since Rh is adjacent to Ru 
in the periodic table but has only one charge state, Rh**. 

From Fig. 6 we conclude that the phonon scattering 
does depend monotonically on P. The maximum value of P 
will be for x=y=0.5 when P,,, , = 2.5 x 10*/cm®. At this value 
Warax = 140. cm K/Watt. This converts to A, =7.1 mW/cm K 
at room temperature. Thus charge transfer scattering can pro- 
duce a very low lattice thermal conductivity in the ruthe- 
nium skutterudites. Whether the carrier mobility and Seebeck 
coefficient can both be kept large at the same time is un- 
known. 

The same sort of thermal conductivity reduction has 
been seen in OsSb,Te and FeSb,Te compounds where one 
finds Os** and Os* of Fe** and Fe**. The scattering strength 
of Fe is comparable to Ru while that of Os is about twice as 
large. These values are, so far based on the analysis of only a 
few samples, and therefore may be modified at a future date. 





Figure 6. 


The Lattice Thermal Resistivity of RhSb, and Various 
Ruthenium-Based Skutterudites at 300K as a Function of the 
Variable-Valence Scattering Parameter, P. 
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Inner Shell Excitations 


The phonon scattering caused by excitation of the elec- 
trons with the f-shell or d-shell of single, isolated transition 
metal atoms has been studied by Slack et al. °*?”. The rare- 
earth ions Sm* and Nd* incorporated into the voids in IrSb, 
exhibit such f-shell phonon scattering™. So far no d-shell 
excitations within single atoms have been seen in 
skutterudites. The f-shell transitions are rather weak phonon 


scatterers compared to the other scattering mechanisms that 
have been discussed previously. 


More Phonon Scattering 
Mechanisms 


Other defects in the skutterudite crystals such as grain 
boundaries, precipitates, and dislocations can produce ther- 
mal conductivity reductions. However we do not believe any 
of these will be useful for controlling the A, in the tempera- 
ture range of interest for thermoelectric devices. 


Conclusions 


Some results on a new family of promising thermo- 
electric materials with the skutterudite crystal structure have 
been presented. These results demonstrate that ZT values of 
the order of unity may be obtainable. Both n-type and p-type 
skutterudites possess good electrical properties although the 
optimum doping levels are much higher for n-type materi- 
als. This is because the electron effective mass is about 12 
times larger than the hole effective mass. The room tempera- 
ture lattice thermal conductivity of the antimony skutterudites 
is between 105 and 160 mW/cm K and needs to be reduced 
by a factor of about 20 in order to achieve high figures of 
merit. Several methods of achieving this have been described. 
The most effective are (a) mass and strain fluctuation; (b) 
void filling; (c) mixed valence electron transfer; and (d) 
charge carrier interactions. The good mechanical properties, 
the low thermal expansion coefficients, and the cubic crystal 
structure are also important if skutterudites are eventually 
employed in good thermoelectric devices. 
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Introduction 


Efficient solid state thermoelectric (TE) devices are 
desirable for several Navy applications, such as shipboard 
air conditioning, cooling of special electronic components 
as well as on-board power generation using thermal sources. 
The performance of a TE device at a temperature T is essen- 
tially related to a thermoelectric figure-of-merit of the mate- 
rial, ZT, used to make the device. The ZT of a thermoelec- 
tric material is described by the relation 


ZT =(o’o/K )T 
where ais the Seebeck coefficient, o is the electrical con- 


ductivity and K is the sum-total of lattice and electronic com- 
ponents of thermal conductivity. The (a0) term is frequently 


referred to as the power factor of the material. A conven- 
tional TE cooling device employs a p-type and an n-type 
segment, (also referred to as legs). Today, a vast majority of 
TE cooling devices employ a p-type leg consisting of a 
Bi, Sb, ,Te, alloy and an n-type leg consisting of a 
Bi,Te, ,.Se, ,, alloy. 

In commercial manufacturing of thermoelectric mod- 
ules, these alloyed materials are prepared as oriented poly- 
crystalline ingots by a Bridgemann technique or a zone lev- 
eling process. The ingots are sliced into appropriate ther- 
moelements, typically, a few mm on each side. To date, the 
best ZT at 300K is about 1.0 for p-type Bi, Sb, Te, alloys 
and about 0.9 for n-type Bi,Te, ..Se, ,, alloys, not significantly 


2.85~ “0.15 
higher than state-of-the-art of thirty years ago [1]. Research 
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on improving the ZT of these bulk thermoelectric materials 
have met with limited success since the 1960's. 

However since 1960's, a variety of thin-film deposi- 
tion techniques such as molecular beam epitaxy (MBE) [2], 
organometallic vapor phase epitaxy (OMVPE) [3] and low- 


pressure chemical vapor deposition (LPCVD) [4] have © 


emerged for depositing a plethora of materials, for applica- 
tions ranging from solid-state lasers to ferroelectrics. These 
thin-film deposition techniques have allowed the synthesis 
of novel materials and enabled material scientists and de- 
vice engineers to apply intriguing concepts such as 
superlattices [5] and quantum-confinement in conventional 
and new devices. However, such techniques of epitaxial 
deposition and the associated analytical techniques for ma- 
terial characterization have rarely been investigated to im- 
prove the performance of thermoelectric materials, until re- 
cently. This is surprising since thermoelectric devices can 
address a substantial commercial and military market, in both 
solid-state cooling and power generation. 


Background 


In a historical note, Whall and Parker [6] suggested a 
possible enhancement in ZT of SiGe-thermoelectric materi- 
als by using superlattice structures, through carrier confine- 
ment in quantum wells to obtain an increased Seebeck coef- 
ficient and by modulation doping to obtain an enhancement 
of carrier mobility. In 1992, a National Thermogenic Work- 
shop was organized [7] to consider various approaches to 
obtain a significant breakthrough in the performance of ther- 
moelectric materials. Thin-film approaches using quantum- 
well structures to enhance power-factor [8, 9] and an argu- 
ment for the use of superlattice structures to decrease ther- 
mal conductivity [10] were proposed in this workshop. The 
motivation for the use of superlattice structures to reduce 
thermal conductivity came from the observed reduction of 
lattice thermal conductivity using solid solution alloying in 
the Bi,Te,- Sb,Te, and Bi,Te,-Bi,Se, material system. In this 
system, the reduction in lattice thermal conductivity has been 
attributed to lattice distortions, from the difference between 
the masses of the constituent atoms, that are effective in 
phonon scattering, thereby impeding the flow of heat. Simi- 
larly, lattice perturbations and the atomic mass fluctuations 
at the superlattice interfaces, are expected to enhance phonon 
scattering rates. The concept of “confined” phonons [1 1] as 
well as a reduction in the thermal conductivity of GaAs/AlAs 
superlattice structures [12], relative to that of AlAs or GaAs 
and depending on the dimensions of the superlattice struc- 
ture, have been experimentally demonstrated. Confined 
phonons imply that they do not propagate suggesting lower 
thermal conduction due to phonons. Also, a planar, thin- 
film, thermoelectric device process technology denoted as 
Monolithically Interconnected, Superlattice-Structured Ther- 
moelements (MISST) [10] was proposed to concurrent-en- 
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gineer the practical application of thin-film, superlattice-type 
structures to high-performance thermoelements. 

Since then, several groups in the United States have 
extended the initial theoretical calculations in Ref. 8 on the 
potential power factor improvement in quantum-well struc- 
tures [13, 14]. Some laboratories have embarked on the ex- 
perimental demonstration of improved ZT at 300K in thin- 
film quantum-well and superlattice structures. 

In this review, we will focus on the experimental re- 
sults in the growth and thermoelectric properties of thin-film 
superlattice and quantum-well structures. A companion pa- 
per by Reinecke et. al in this issue discusses the theory be- 
hind some of the thermoelectric properties of superlattices. 
We will first describe some aspects of thin-film deposition 
of superlattice structures in thermoelectric materials. Fol- 
lowing this, we describe some of the key results that have 
recently emerged with the use of quantum-confined and 
superlattice structures. Various approaches to the fabrica- 
tion of a thin-film TE device will be described briefly. Fi- 
nally, we discuss the major research issues yet to be fully 
addressed toward the eventual implementation of thin-film 
thermoelectric materials in practical, solid-state cooling and 
power generation systems. 


Deposition of 
Thermoelectric 
Superlattice Structures 


In this section we will consider two studies of epitaxial 
growth and characterization of thermoelectric properties of 
superlattice structures, where such properties have been in- 
vestigated (and reported) as a function of superlattice dimen- 
sion. Such a dimensional-related quantification is impor- 
tant to distinguish the improved properties from a superlattice 
approach, rather than arising purely from an epitaxial depo- 
sition process. 

One study is on n-type, MBE-grown PbTe/Po, Eu Te 
and related superlattice structures, grown on semi-insulat- 
ing BaF, substrates. This system has been originally studied 
and reported by Harman et. al of MIT Lincoln Laboratories 
[15]. Second, we will describe our work on thermoelectric 
properties of p-type, OMVPE-grown Bi,Te,/Sb,Te, 
superlattice structures on semi-insulating GaAs and Al1,O, 
substrates [16]. Farmer et. al of Lawrence Livermore Na- 
tional Laboratory have reported on sputter-deposited 
Bi, Sb, /X (X=Bi,Te,, Bi and PbTe, Se, .) superlattice struc- 
tures on semi-insulating Al,O, substrates [17]. However, their 
early results on power-factor of these sputter-deposited 
superlattice thin-films are inferior or about comparable to 
non-superlattice films of pure Bi and Bi, Sb, ,. Elsner et. al 
of Hi-Z Technology, in cooperation with Naval Research 
Laboratory, have reported on some significant enhancement 





of thermoelectric properties of Si/Si, Ge, , superlattice struc- 
tures grown on conducting Si substrates [18]. The conduct- 
ing nature of the Si substrate, especially when it is lightly- 
doped (therefore has a high Seebeck coefficient) and is of 
the same polarity as the epitaxial film that is being investi- 
gated, makes the interpretation of properties of superlattices 
very difficult. 


MBE of PbTe / Pb, Eu Te 
Superlattices 


The MBE growth of PbTe and its alloys, including Pb, 
,Eu,Te and Pb, Eu Te Se, ., have been researched in the past 
[19] for use in lead-salt diode lasers and detectors. Their 
growth have been demonstrated on single-crystal, cleaved 
BaF, substrates. The insulating nature of the BaF, substrates 
is very convenient for studying the thermoelectric properties 
as well. Harman et al. [15, 20] have pursued this materials 
system to evaluate the predicted quantum-well effects in 
PbTe/Pb, Eu Te system and its advantages for ZT improve- 
ment. 

Typical PbTe/Pb, Eu Te (x~0.073) multi-quantum well 
structures are grown on BaF, substrates with a 2000 A Pb Eu,_ 

yle buffer. PbTe is the quantum- -well material, typically 20 
to 50 A thick. Pb, Eu Te, which is the barrier to provide the 
two-dimensional quantum- -confinement of electrons in the 
PbTe well, is about 460 A thick [15, 20]. The PbTe layers 
are left undoped and the barriers are doped with Bi [15]. 
This approach to doping is called modulation doping [21]. 
More details on the growth and charactetrization of these 
structures will be published soon [22]. 


OMVPE of Bi,Te, / Sb,Te, 
Superlattices 


To our knowledge, there had been no reported OMVPE 
growth of Bi,Te, and related materials about two years ago 
when we began the experimental evaluation of thermoelec- 
tric properties of Bi,Te,/Sb,Te, superlattices. However, 
OMVPE, being a kinetically- limited deposition process rather 
than an equilibrium-driven process, has been successful in 
the synthesis of a variety of materials including those with 
miscibility gaps. Further, the availability of a wide range of 
organometallic precursors for Te with extremely low crack- 
ing temperatures was a good starting point. 

A schematic of the thin-film Bi,Te /Sb,Te, superlattice 
structure that we have proposed as an alternative to the con- 
ventional bulk alloy is shown in Figure 1. In this review, we 
will consider the thermoelectric properties of two types of 
superlattice structures. One type is the symmetric 
superlattice, where the thickness of the Bi,Te, and the Sb,Te, 
layers are kept the same. The other is the non-symmetric 
superlattice, where the thickness of the Sb,Te, layer is larger 


than the Bi,Te, layer. 

OMVPE/LPCVD growth system used to deposit the 
materials, has been designed to allow the fabrication of high 
quality superlattice structures. The OMVPE growth process 
is easily amenable to scaling for large-area applications, if a 
thin-film device fabrication demands it, as discussed in the 
following section on approaches to thin-film thermoelectric 
devices. We have grown Bi,Te, and Sb,Te, films on both 
semi-insulating GaAs substrates and Al 0, substrates [23]. 
The growth has been achieved by a novel combination of 
high-temperature cracking of organometallic precursors and 
low substrate temperatures [24]. This new approach to 
OMVPE has been instrumental in obtaining single-crystal- 
line Bi,Te, films and growth of thermally stable, single-crys- 
talline films of Sb,Te,. 

X-ray diffraction data obtained on Bi,Te, films on Al,O, 
and GaAs indicate single-crystallinity. Low- -energy electron 
diffraction data has confirmed the single-crystallinity of the 
film as well. The x-ray full-width at half-maximum, a qual- 
ity parameter, of the films are comparable with a recently 
reported MBE-grown Bi,Te, film on a Sb,Te, substrate [25], 
in spite of a larger lattice mismatch between the Bi,Te, film 
and an ALO, substrate visa-vis Bi,Te, film and a Sb ‘Te, sub- 
strate. The ‘stoichiometry of Bi, Te, ‘and Sb,Te, films have 
been confirmed by Rutherford backscattering ‘and X-ray pho- 
toemission spectroscopy (XPS), respectively. The as-grown 
Sb,Te, films are always p-type, determined by Hall-effect 
measurements. The as-grown Bi,Te, films are always n-type, 
while the Bi,Te/Sb,Te, superlattice structures were always 
found to be p-type. ‘More details of the growth of these ma- 
terials are presented elsewhere [26]. 

An important concern in the growth of superlattices, 
in any material system, is the possibility of inter-mixing be- 
tween the constituent layers of the superlattice. In this re- 





Figure 1. 


Schematic of a short-period p-type Bi,Te /Sb, Te, superlattice 
structure that we have investigated as an alternative to a 
conventional bulk alloy in p-type (Bi, ,,Sb, ,,),Te,, In this paper, 
when L.=L,, it is denoted as a symmetrical superlattice and 
when L.>L,, it is denoted as a non-symmetrical superlattice. 
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gard, it is instructive to examine the crystal structure of a 
layered Bi,Te,/Sb,Te, material, along the growth direction. 
The schematic of this crystal structure for a~17 A superlattice 
dimension is shown in Figure 2. We note the periodic van 
der Walls-like bond within Bi,Te, and between Bi,Te, and 
Sb,Te, layers in this structure. 





Figure 2. 


Crystal structure of a layered Bi,Te/Sb,Te, material along 
the growth direction. Note the presence of van der Walls-like 
bonding along the c-axis. 
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It has been predicted that this type of layered struc- 
ture, with van der Walls bonding along the growth direction, 
lends itself amenable for obtaining abrupt surface termina- 
tions [27] thereby potentially leading to high-quality 
superlattice interfaces. This was evident in the XPS spectra 
obtained on a very thin (< 50 A) Sb, Te, layer grown on a 
thick Bi,Te, layer. XPS data, evaluated in both normal and 
glancing modes, indicated a negligible presence of Bi in the 
surface Sb,Te, layer. This provided some evidence that abrupt 
termination of the Bi,Te, layer, prior to the growth of Sb,Te, 
layer, was likely achieved. The quality of the Bi,Te JSb, Te, 
superlattice structures were assessed by single- crystal x-ray 
diffraction. In the diffraction spectrum, clear satellite peaks 
(corresponding to the periodicity of the superlattices) [28] 
were obtained. Good agreement was obtained between the 
nominal value of the period estimated from growth-rates of 
the individual layers and the superlattice period deduced from 
the position of the satellite peaks. We have also recently evalu- 
ated these superlattices by double-crystal x-ray diffraction 
and transmission electron microscopy (TEM) [26]. The TEM 
data, considering that the sample preparation for microscopy 
with the “soft” Bi,Te, and related materials are in its infancy, 
has indicated that good quality interfaces were achieved in 
the OMVPE-grown Bi,Te,/Sb,Te, superlattices. 
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Thermoelectric 
Properties of Quantum- 
well and Superlattice 
Structures 


Power Factor Enhancement in 


PbTe / Pb, ,,-Fu, ,,,re Quantum 
Wells 


Harman et al. have recently reported on the in-plane 
Seebeck coefficient in PbTe quantum wells, as a function of 
free-carrier concentration, for a pair of quantum-well dimen- 
sions [15] and compared them to bulk and epitaxial PbTe 
materials. This data is reproduced in Figure 3. The excel- 
lent agreement of Seebeck coefficient as a function of dop- 
ing level in epitaxial materials, visa-vis those of bulk mate- 
rials, attests to the validity of the measurement of Seebeck 
coefficients in thin epitaxial films. The Seebeck coefficient 
of quantum-well structures show a significant enhancement 
over bulk materials only when quantum-well dimensions are 
about 20A as opposed to 40 A. This experimental investiga- 
tion was done to test the theoretical model [8, 20], which 
predicts a large increase in power factor (and hence ZT) in 
multiple quantum-well structures. Good agreement was 
found between the experimental results and the calculations 





Figure 3. 


Seebeck coefficient of bulk PbTe and PbTe quantum-wells 
as a function of carrier concentration [15]. 
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[20], giving an apparent experimental confirmation of the 
model. The ability to obtain higher Seebeck coefficient, for 
the same carrier density with apparently no loss of electronic 
mobility, translates to improved power factor for the 20-A 
quantum-well structure. The power factor of the 20-A PbTe 
quantum-well is almost a factor three better than that of bulk 
PbTe. Some recent work indicates that improved growth pa- 
rameters can lead to higher Seebeck coefficient (also) for 
the 40-A quantum-well compared to bulk PbTe. 


Potential ZT in PbTe/ 
Pb Eu u, o734& Quantum-Well 
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The improved power factor in the 20-A PbTe quan- 
tum-well, in combination with an estimated thermal conduc- 
tivity, results in a Z,,T of 1.3 at 300K. This is a factor of 
three higher than that of the best bulk PbTe materials [15]. 
The Z,,,7 is a two-dimensional figure-of-merit for the two- 
dimensional PbTe quantum-well and neglects the parasitic 
thermal conductivity of the barrier layer; when this is ac- 
counted for, a three-dimensional Z,,T of 0.11 is estimated 
[15, 29]. 


Power Factor Enhancement in 
Bi,Te, / Sb,Te, Superlattices 


Hall-data indicates the residual hole concentration in- 
creases in the symmetrical Bi,Te,/Sb,Te, superlattice struc- 
tures, with smaller superlattice dimension. Note here that 
the superlattice dimension refers to the thickness of each of 
the Bi,Te, and Sb,Te, components, which is kept about equal 
for the symmetrical ‘superlattice structures. The noticeable 
sharp increase in residual p-type behavior with smaller 
superlattice dimension (achieved by more Bi,Te, and Sb,Te, 
growth switchings) suggests that Te-vacancy defects are pos- 
sibly formed at the superlattice interfaces. Such defects are 
expected to behave as acceptors [30]. Recently, we have 
minimized this rapid increase in carrier concentration in short- 
period superlattices with higher Te-vapor pressure during 
OMVPE growth, consistent with this Te-vacancy argument. 

The behavior of the Hall hole-mobility as a function 
of the superlattice dimension, for the case of symmetrical 
superlattice structures, is indicated in Figure 4a. Here, we 
notice that for all the superlattice dimensions, the hole mo- 
bility is higher than that of a corresponding (Bi, Sb, ,),Te, 
alloy with a similar resistivity. Even for a superlattice di- 
mension of ~10A, the mobility of the superlattice structure 
is better than the comparable solid solution alloy. This is 
evidence of avoiding or minimizing alloy scattering of car- 
riers in Bi,Te,/Sb,Te, superlattices, relative to the 
(Bi, ,Sb, ,), Te, alloy. However, we do notice a big drop in 


mobility when the superlattice dimension falls below the unit 
cell dimension of ~30 A. The dependence of carrier mobil- 
ity on the superlattice dimension for the non-symmetrical 
superlattice structures, with a constant Bi,Te, layer thick- 
ness of ~15 A, is shown in Figure 4b. We. once again note 
that for a Sb,Te, layer dimension >30 A, the mobility of the 
superlattice structure exceeds that of a typical solid solution 
(Bi, Sb, .) Te, alloy. 

The dependence of the 300K Seebeck coefficient of 
the symmetrical and non-symmetrical Bi,Te,/Sb,Te, 
superlattice structures on the superlattice-dimension is shown 
in Figure 5a and 5b, respectively. For the symmetrical case, 
we see a near 50% increase in the Seebeck coefficient with 
nearly a ten-fold reduction in the superlattice dimension. This 
increase is relatively small compared to a near 100% increase 
in Seebeck coefficient, in going from a 40-A well to a 20-A 
well, in the case of PbTe/Pb, ,,,Eu, ,,, Te quantum-well sys- 
tem (Figure 3). 





Figure 4. 


Carrier mobility variation with superlattice dimension in (a) 
symmetrical and (b) non-symmetrical Bi,Te/Sb,Te, 
superlattice structures. 
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Figure 5. 


Variation of the Seebeck coefficient with superlattice 
dimension in (a) symmetrical and (b) non-symmetrical 
superlattice structures. 
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The above difference in the behavior of Seebeck coef- 
ficients with the superlattice dimension can be interpreted in 
terms of the likely difference in the properties of Bi,Te,/ 
Sb,Te, and PbTe/Pb, ,,,Eu, ,,,Te hetero-interfaces. The 
bandgap difference between PbTe and Pb, ,.,Eu, ,,, Te is about 
0.31 eV at 300K, leading to a conduction-band offset of ~0.17 
eV in the n-type PbTe/Pb, ,,,Eu, ,,,Te hetero-junction [20]. 
This conduction-band offset provides the necessary barrier 
for the quantum-confinement of electrons in the PbTe-well. 
In contrast, the bandgap difference between Bi,Te, and Sb,Te, 
materials is expected to be less than 0.05 eV at 300K [10]. 
The distribution of this small bandgap difference in the Bi,Te./ 
Sb,Te, superlattice system, between a valence-band offset 
and a conduction-band offset, has not been measured. How- 
ever, if we assume that these band-offsets are about equal 
(~0.03 eV), we expect only a weak hole-quantum-confine- 
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ment in the p-type Bi,Te /Sb,Te, superlattices at 300K. This 
could explain the slow rise in Seebeck coefficient with smaller 
superlattice dimensions in the Bi,Te,/Sb,Te, system. We also 
note a small drop in Seebeck coefficient with decrease in 
superlattice dimension, in the range of 20 A to ~10 A, prob- 
ably from a result of the sharp increase in carrier levels. The 
300K Seebeck coefficients of the non-symmetrical Bi,Te,/ 
Sb,Te, superlattice structures also increase slowly with a re- 
duction in the superlattice dimension as seen in Fig. 5b. 

It is important to note in both Fig. 5a and 5b, the 
Seebeck coefficients of the short-period symmetrical and non- 
symmetrical superlattice structures are higher than that of a 
BiSbTe,-alloy film and a typical (Bi, .Sb, ,),Te,-alloy film, 
respectively. This can be attributed in part to the Fermi-level 
differences between the respective superlattice structure and 
the corresponding alloy film. The relation to obtain the 
Seebeck coefficient of a p-type material is given by the rela- 
tion below: 


jf Eo(E) dE 
k E=E 


a= 


Here, © is the conductivity, E* is the reduced energy 
of holes in the valence band and E,* is the reduced Fermi- 
energy. For the same conductivity, the higher mobility in 
the superlattice (compared to the alloy) leads to a lower car- 
rier concentration and therefore a lower hole-E,*. This can 
potentially lead to a higher Seebeck coefficient for the 
superlattice structure when compared to the alloy. 

The effect of an improved hole mobility for a larger 
Seebeck coefficient is to increase the power factor in short- 
period p-type Bi,Te,/Sb,Te, superlattice structures. We have 
observed that that the power factor of the non-symmetric 
superlattice structures are consistently higher than the sym- 
metrical superlattices, due to the effect of higher mobilities 
for the same resistivities. We have observed a best power 
factor in the range of 55 uW/cm-K? in a non-symmetrical 
Bi,Te,/Sb,Te, (15 A/45 A) superlattice structure at 300K. 

An important concern in the measurement of thermo- 
electric properties of thin-film is the possible effect of the 
substrate. We have compared the Hall-effect carrier concen- 
tration and mobility on Bi,Te,/Sb,Te, symmetrical 
superlattices (for a constant superlattice dimension) on a 
semi-insulating GaAs and an insulating ALO, substrate. The 
comparable data on the two substrates, given our observa- 
tion that the morphological-quality of Bi,Te, and Sb,Te, films 
are generally better on GaAs substrates than on Al,O, sub- 
strates, suggests to the absence of any effects of the sub- 
strate. Also in this regard, we have grown several thin-film 
p-type (Bi, ,.Sb, ,,), Te, alloy layers on semi-insulating GaAs 
substrates and compared their Seebeck coefficients to bulk 
p-type (Bi, ,.Sb, ,.),Te, alloy materials, as a function of free- 





carrier levels. We have found good correlation [16] of the 
Seebeck coefficient, as a function of the free-carrier level, 
between these epitaxial thin-films and those of bulk materi- 
als [30]. This further attests to the validity of our measure- 
ment of the Seebeck coefficient in thin-films, similar to the 
above-discussed correlation in Seebeck coefficients of epi- 
taxial-PbTe and bulk-PbTe materials. 


Thermal Conductivity Reduction 
in Bi,Te,/Sb,Te, Superlattices 


An important objective of our experimental research 
effort on short-period Bi,Te,/Sb,Te, superlattices has been 
to demonstrate that phonon-scattering at the superlattice in- 
terfaces can reduce the lattice thermal conductivity and then 
verify if it is indeed possible to obtain values lower than in 
solid-solution (Bi, ,.Sb,,,),Te, alloys [10] from additional 
phonon-scattering processes unique to periodic structures. 
For example, in a superlattice, the minimum phonon energy 
required to produce an Umklapp process is reduced relative 
to a bulk alloy. Thus a measurement of thermal conductivity 
is important to take further advantage of the improved elec- 
trical transport and power factor of Bi,Te,/Sb,Te, 
superlattices, in the overall optimization of ZT values. 

The modified 3-@ technique [31] has been used to 
measure thermal conductivity (A) of thin-films from an indi- 
rect measurement of thermal diffusivity. Essentially, a thin 
metal strip with four pads to measure current and voltage is 
deposited on the film surface. The metal strip serves as the 
heater (by I? R heating) and a thermometer because of the 
temperature-dependent electrical resistance of the strip is 
related to the thermal conductivity of the thermoelectric 
material. Thus the thermal conductivity can be related to 
measured voltages. The application of an ac current (I) with 
an angular frequency @ causes a I’ R heating temperature 
wave of frequency 2@ to diffuse into the substrate. Since the 
film resistance itself varies with temperature, the resistance 
oscillation at 2 multiplied by the excitation current at @ 
produces a voltage oscillation at 3m. The 3@ voltage can be 
conveniently measured by a lock-in amplifier technique. 

Figure 6 shows the thermal conductivity data on both 
symmetrical and non-symmetrical Bi, Te,/Sb,Te, superlattice 
structures as a function of the superlattice dimension. The 
thermal conductivity values apparently come down with a 
shorter superlattice dimension, although the superlattice 
structures in general (for periods in the range of ~100 A or 
less) do indicate reduced thermal conductivities compared 
to bulk materials. We observe an apparent minimum thermal 
conductivity when the Bi,Te, or Sb,Te, dimension approaches 
~30 A, the unit cell dimension. For these superlattice di- 
mensions, we obtain thermal conductivity values as low as 
~2.0 mW/cm-K. It is worth observing that both the sym- 
metrical and non-symmetrical superlattices apparently of- 
fer about the same value for minimum thermal conductivity 


at a dimension of ~30A, indicating that the thermal conduc- 
tivities in these superlattices are likely dominated by phonon 
scattering and reduced phonon group velocities as a result 
of the superlattice structure, rather than the bulk-like phonon 
scattering processes. It is also interesting to note that when 
the superlattice dimension is less than ~30 A, the unit cell 
dimension, the measured thermal conductivities actually in- 
crease and tend to approach the values measured for the re- 
spective alloy films. 





Figure 6. 


Variation of thermal conductivity with superlattice dimension 
in symmetrical and non-symmetrical Bi, Te /Sb, Te, superlattice 
structures. 
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Comments on 3w-Measured 
Thermal Conductivities in 
Superlattices 


It is worth noting that the measured 3@ thermal con- 
ductivity value for the p-(Bi, ,.Sb, ,.),Te, alloy film is ~16.5 
mW/cm-K and for the p-(BiSb) Te, alloy film is ~10.7 
mW/cm-K. These are approximately equal to the geometric 
mean of known thermal conductivities, along the “a” and 
“c” directions, in bulk crystals of these materials [1, 30]. 
Thus the 3@-thermal conductivity method, at least for the 
frequencies (~1000 Hz) used in this study, could also be af- 
fected by the spread of the thermal wave in the direction 
parallel to the superlattice interfaces. It appears that the 3 
thermal conductivity value is not likely to be dominated solely 
by the component perpendicular to the superlattice interfaces. 
This was further evident in the observation that the 3@-mea- 
sured thermal conductivities in GaAs/Al, Ga, ,As 
superlattices were comparable to the in-plane thermal con- 
ductivity values of GaAs/AIAs superlattices measured by AC 
calorimetry [12] and larger than the perpendicular thermal 
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conductivity values in GaAs/AIAs superlattices, determined 
by a pico-second reflectivity-optics method [32]. For ex- 
ample, the pico-second optics method on a ~50A/50A GaAs/ 
AlAs superlattice indicated a factor-of-2 to 3 lower thermal 
conductivity compared to values obtained by the 3@ method 
or the in-plane method for similar superlattice dimensions. 

In essence, we note that the thermal conductivity val- 
ues of the Bi Te /Sb, Te, superlattice structures need to be 
evaluated by a “true’ * through- -thickness method (such as in 
Ref. 32) to measure the component perpendicular to the in- 
terfaces and a “true” in-plane method, with the substrate re- 
moved as in Ref. 12, to obtain the thermal conductivity par- 
allel to the superlattice interfaces. 


Potential ZT in Bi,Te,/Sb,Te, 
Suyerlattice Structures 


The power factors in Bi,Te,/Sb,Te, superlattice struc- 
tures, in this work, have been measured i in a direction paral- 
lel to the superlattice interfaces. As mentioned in the previ- 
ous section, the directional nature of the 3-@ thermal con- 
ductivity value is not completely clear at this point in time. 
In addition, unlike GaAs and AlAs materials, Bi,Te, and 
Sb,Te, materials do exhibit anisotropic thermal conductivi- 
ties. This makes the determination of ZT values of the Bi,Te/ 
Sb,Te, superlattice structures not straightforward. 

However, as noted earlier, in the short-period Bi,Te,/ 
Sb,Te, superlattice structures where the thermal conductiv- 
ity values are as low as ~2.0 mW/cmK, the superlattice in- 
terfaces rather than the bulk properties probably control the 
thermal conductivity values. So, we need to potentially cor- 
rect for typical anisotropy of thermal conductivities in 
superlattice structures, rather than for bulk-like anisotropies. 
We can estimate the possible anisotropy of thermal conduc- 
tivity values in Bi,Te,/Sb,Te, superlattice structures based 
on some recent comparative measurement [33] of thermal 
conductivities in GaAs/AIAs superlattices, along both per- 
pendicular and parallel to the superlattice interfaces. This 
work in Ref. 33 suggests that, in a worst-case scenario, the 
thermal conductivity value parallel to superlattice interfaces 
could be higher than that along perpendicular to the inter- 
faces by a factor of about 2.6; similarly, we estimate that the 
thermal conductivities of the Bi, Te /Sb,Te, superlattices could 
be probably higher in the plane of the superiattice- interfaces 
by this ratio. This leads us to estimate a three-dimensional 
ZT,,, of around 2.0 at 300K for some of the best Bi,Te,/Sb,Te, 
superlattice samples in the direction parallel to the superlattice 
interfaces. This data is indicated in Table 1. In this estima- 
tion, we have used the power factors measured at MIT-Lin- 
coln Laboratories for the Bi,Te,/Sb,Te, superlattices. Also 
shown in Table | are the ZT values in comparable alloy films, 
based on RTI-measured power factors. We note that the 
power factors measured at RTI on the Bi,Te,/Sb,Te, 
superlattices, done similar to these alloy samples, have agreed 
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well with the measurement at MIT-Lincoln Laboratories. 
Hence, the ZT enhancement in Bi,Te,/Sb,Te, superlattice 
structures over comparable alloy films, grown by OMVPE, 
and the state-of-the-art bulk materials (typically around 0.9 
to 1.0) are noteworthy. 

The estimation of ZT,, for the Bi,Te,/Sb,Te, 
superlattices, in the direction perpendicular to the superlattice 
interfaces, is much more complicated and will require the 
measurement of Seebeck coefficient and electrical conduc- 
tivity in this direction. Such an effort is in progress. 


Approaches to Thin-Film 
Thermoelectric Devices 


Recently, Farmer et. al [17] have denionstrated a min- 
iature thermoelectric cooling device with an 8-micron-thick, 
sputtered, p-(Bi,,.S “ Rah Te, and n- Bi, ,Sb, , thin-films, re- 
spectively, on 25-um- -thick mica substrate. "These devices 
were operated with a duty-cycle of 0.3, instead of a conven- 
tional, continuous, direct-current operation and a cooling of 
5 to 10°C below the ambient was obtained. This is an en- 
couraging result and portends well for thermoelectric devices 
based on high-ZT superlattice thin-films. 

The results on superlattice and quantum-well struc- 
tures grown by epitaxial techniques for high-ZT thermoelec- 
tric applications, although appear very promising, need to 
be balanced with the fact that these epitaxial techniques are 
ideally suited for producing thin layers and structures, typi- 
cally ~ 5 to 10 um thick, as a result of smaller growth rates. 
Typically, growth rates are 0.1 m/min in OMVPE systems 
and 0.02 m/min in MBE systems. These growth rates are 
adequate for the deposition of thin-films required for many 
electronic applications. However, practical thermoelectric 
elements require 200- to 1000-{um-thick layers [1] in the di- 
rection of heat flow, i.e., along the thermal gradient. The 
Peltier contact area is used to pump heat from the cold junc- 
tion to the hot junction. Parasitic thermal conduction from 
the hot junction toward the cold junction and Joule heating 
oppose the Peltier-induced heat flow. A smaller thermal gra- 
dient is required to reduce this reverse flow of heat i.e., main- 
tain certain optimum spatial separation between the hot and 
cold junctions. Typically, this separation is ~ 200 to 1000 
tum and it would seem that it will be difficult to grow high- 
quality epitaxial superlattice structures to such thickness, 
even by present-day thin-film deposition techniques like 
MBE and OMVPE. 

One approach [34] to utilizing the superlattice struc- 
tures as high-performance thermoelectric elements, there- 
fore, may be to process the thin-films (~5 to 10 pum) into a 
thermoelectric element of effective lengths of ~100 to ~500 
lum using a planar monolithically-interconnected, 
superlattice-structured thermoelement (MISST) approach 
[10]. With the MISST approach, the ohmic contacts can po- 





tentially generate significant I’R heating, leading to a degra- 
dation of performance from that obtainable with an enhanced 
ZT. Initial modeling of the cooling-performance of the 
MISST device indicates that using only four segments on 
each of the n- and p-type element, with an epitaxial layer 
thickness of ~10 um, a resistivity of ~3E-3 ohm-cm for the 
superlattice material, a Peltier contact resistivity of 1E-6 
ohm-cm’, an ohmic-contact resistivity of SE-7 ohm-cm’, an 
effective device ZT of ~1.5 may be obtained with an intrin- 
sic ZT of ~2.3. Therefore, good quality ohmic contacts at 
each of the vertical segment are necessary to reduce  R 
losses, which can affect the overall cooling efficiency. We 
note that the use of a limited number of segments in the 
MISST device in conjunction with a reduction in radiative 
heat loss from the ohmic-contact surfaces, and minimiza- 
tion of thermal conduction losses through the supporting- 
superstrate like Kapton, would help ameliorate the parasitic 
thermal exchange with the ambient. 

The MISST concept, if successful, can alleviate the 
necessity of epitaxially growing thick epitaxial superlattice- 
type layers, as required in other possible approaches to fab- 
rication of a thin-film thermoelectric device [9]. The MISST 
device processing technique can enable the cost-effective 
fabrication of high-performance thermoelements from 
superlattice-type materials as well as non-superlattice thin- 
film materials. Fairly simple, conventional microelectronic 
processing is required to fabricate these thermoelements. Two 
sets of wafers need to be processed, one for the n-type leg 


and the other for a p-type leg, which are then interconnected 
to form a cooling element. The attachment of the thermally 
insulating superstrate like Kapton will enable the removal of 
the growth-substrate. This superstrate provides mechanical 
support for the interdigitated n-type and p-type legs. While 
a smaller thickness of the Kapton superstrate is preferable 
for reducing the parasitic thermal conduction loss, a larger 
thickness will facilitate device handling. 


Near-Term Research 
Issues 


In earlier sections, we pointed out several aspects of 
thermoelectric properties of both PbTe quantum-well struc- 
tures and Bi,Te,/Sb,Te, superlattices that merit a closer look. 
Specifically, in MBE-grown PbTe/Pb, ,,.Eu, ,,Te structures, 
the requirement of a fairly thick (460 A) Pb, ,,,Eu, ,,, Te bar- 
rier compared to the quantum-well (20 A) results in a sub- 
stantial degradation of an estimated, effective 3-dimensional 
ZT compared to an excellent ZT,,, of ~1.2 for the PbTe quan- 
tum-well. There is a clear need to measure thermal conduc- 
tivity of these PbTe/Pb, ,,.Eu, ,,,Te structures, by methods 
such as 3-@, to obtain an estimate of Z,,T. The thermal con- 
ductivity results on Bi,Te,/Sb,Te, superlattice structures [16] 
should be an incentive to explore these PbTe/Pb, ,,.Eu, ,,, Te 
structures as well. As suggested by Harman et. al [15], other 
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quantum-well/barrier combinations such as Bi Sb, /Pb Eu, 
, Te Se, , or a suitable alternative need to be explored to ob- 
tain higher Z,.T of up to ~8.0 and correspondingly Z,,T of 
near 3.0. This should help overcome some of the restric- 
tions posed by the requirement of thick barrier layers in the 
quantum-well approach. 

The early results on thermoelectric properties of Bi,Te / 
Sb,Te, superlattice structures appear to suggest that it may 
be possible to break the traditional, conflicting requirements 
for high electrical conductivity and low thermal conductiv- 
ity encountered in bulk thermoelectric materials. 
Superlattices may offer a route to minimize carrier scatter- 
ing and obtain improved transport properties including ther- 
moelectric power-factor, by avoiding solid-solution alloy- 
ing; yet the superlattice interfaces can enhance phonon scat- 
tering rates and reduce lattice thermal conductivity. One of 
the key issues in this material system remains - what are the 
conduction and valence band offsets so that a better theoreti- 
cal understanding of the ot served enhancement in thermo- 
electric properties can emerge. Another key issue is the ac- 
curate measurement of the thermal conductivity in 
superlattice structures, along specific directions. 

From an experimental standpoint, it is worth noting 
the importance of high-quality Bi,Te,/Sb,Te, superlattice 
structures on improvement in thermoelectric properties. We 
have observed that the effect of an intentional disorder in the 
superlattices is to reduce the Seebeck coefficient and increase 
the thermal conductivity. The disorder was simulated by in- 
tentionally making the Bi,Te, component of the superlattice 
as an alloy layer. Such controlled studies of the Bi,Te /Sb,Te, 
superlattice in conjunction with measurements of hetero-in- 
terface properties are urgently needed. A maior issue re- 
mains the reliable measurement of all three thermoelectric 
properties of thin-film structures, namely the Seebeck coef- 
ficient, the electrical conductivity, and the thermal conduc- 
tivity, both parallel and perpendicular to the superlattice in- 
terfaces. Uher and Moreli discuss some of these aspects in a 
companion paper in this issue. 

Development and testing of thin-film thermoelectric 
devices employing some of the superlattice structures (with 
the higher estimated ZT values) are needed. Thin-film de- 
vice design optimization will include identification of the 
necessary dimensions based on thermoelectric material pa- 
rameters. Significant process development issues will have 
to be addressed depending on the thin-film device approach 
employed in the fabrication. The device process develop- 
ment such as joining the n-type and p-type superlattice ele- 
ments together, aligning and joining them to an electrical 
path, and providing for the structural rigidity of the resulting 
structure, are expected to be significant engineering prob- 
lems that can be overcome once we are convinced that we 
do have ZT>1 in thermoelectric materials based on 
superlattice structures. 


40 Naval Research Reviews 


Acknowledgments 


The author acknowledges the technical contributions 
of Mr. Thomas Colpitts, Mr. Eric Watko, Mr. Mike Mantini 
and Mr. Dave Malta during the course of this work. The 
author also acknowledges Dr. Ted Harman of MIT Lincoln 
Laboratories for the verification of power factors in the 
Bi,Te,/Sb,Te, superlattice structures as well as many valu- 
able comments during the review of this manuscript. 


The editor wishes to acknowledge Dr. Ted Harman of 
MIT Lincoln Laboratory for assistance in the review of the 
scientific and technical details of the article and for volun- 
teering to assist in measuring properties of thin film struc- 
tures reported in this article. 


Biography 


Dr. Rama Venkatasubramanian received his Ph.D. in 
Electrical Engineering in 1988 from Rensselaer Polytechnic 
Institute, Troy, New York. He was a recipient of the Allen B. 
Dumont Prize at Rensselaer for outstanding academic 
achievement. Prior to his graduate work at Rensselaer, he 
obtained his B.S in Electrical Engineering in 1983 from the 
Indian Institute of Technology, Madras, India. He is a re- 
cipient of the National Talent Search Scholarship during his 
undergraduate work. Since 1989, he has been with the Cen- 
ter for Semiconductor Research, Research Triangle Institute, 
in North Carolina. Presently, he is a Senior Research Engi- 
neer working in the area of thermoelectrics, photovoltaics 
and other electronic materials for device applications. He 
has over 75 publications in the area of heteroepitaxy, char- 
acterization of electronic materials, novel optoelectronic 
materials based on quantum-confined structures, 
thermoelectrics, photovoltaics and other devices. 


REFERENCES 


FD. Rosi and E.G. Ramberg, “Evaluation and 
Properties of Materials for Thermoelectric Applica- 
tions,” in Thermoelectricity, Ed. by P. H. Egli, John 
Wiley and Sons, NY (1960). 

A.Y. Cho and J.R. Arthur, Progress in Solid State 
Chemistry, G. Somorjai and J. McCaldin Eds. 10, 
(1973), p.157. 

G.B. Stringfellow, Organometallic Vapor Phase 
Epitaxy: Theory and Practice, Academic Press, 
London, (1989). 

M.E. Jones and D.W. Shaw, Growth from the Vapor, 
from Treatise on Solid State Chemistry, 5, N.B. 
Hannay, Ed., Plenum, New York (1976). 

L. Esaki and R. Tsu, IBM J. Res. Dev., 14, 61 





(1970). 

T.E. Whall and E.H. C. Parker, Proc. of First Euro- 
pean Conf. on Thermoelectrics, D.M. Rowe, Ed., 51 
(1987) and references there in. 

lst National Thermogenic Cooler Conference, Sept. 
17, 1992, Fort Belvoir, VA. 

L.D. Hicks and M.S. Dresselhaus, “The Effect of 
Quantum Well Structures on the Thermoelectric 
Figure of Merit,” Presented at the Ist National 
Thermogenic Cooler Workshop, Sept. 17, 1992, Fort 
Belvoir, VA. 

T.C. Harman, “PbTeSe/BiSb Short-Period 
Superlattice as a New Thermoelectric Cooler 
Material,” Presented at the 1st National Thermogenic 
Cooler Workshop, Sept. 17, 1992, Fort Belvoir, VA. 
R. Venkatasubramanian, “Organometallic Epitaxy of 
Bi,Te, and Related Materials and the Development of 
Planar, Monolithically-Interconnected, Superlattice- 
Structured, High-Efficiency Thermoelectric Ele- 
ments,” Presented at the 1st National Thermogenic 
Coler Workshop, Sept. 17, 1992, Fort Belvoir, VA. 
B. Jusserand and D. Paquet, “Confined and 
Propogative Vibrations in Superlattices,” in 
Heterojunctions and Semiconductor Superlattices, 
Ed. by G. Allan, G. Bastard, M. Launoo, and M. 
Voos, Springer Verlag, 1986. 

T. Yao, Appl. Phys. Lett., 51, 1798 (1987). 

D.A. Broido and T.L. Reinecke, Appl. Phys. Lett., 
67, 1170 (1995). 

J.O. Sofo and G.D. Mahan, Proc. of the 13th Interna- 
tional Conf. on Thermoelectrics, AIP Conf. Proc., 
316, 239 (1995). 

T. Harman, Proc. of ARPA Thermoelectric Materials 
Workshop, Arlington, VA, November 1995. 

R. Venkatasubramanian, Proc. of ARPA Thermoelec- 
tric Materials Workshop, Arlington, VA, November 
1995. 

J. C. Farmer, Proc. of ARPA Thermoelectric Materi- 
als Workshop, Arlington, VA, November 1995. 

N.B. Elsner, S. Ghamaty, J.H. Norman, J.C. Farmer, 
R.J. Foreman, L.J. Summers, M.L. Olsen, P.E. 
Thompson, and K. Wang, Proc. of the 13th Interna- 
tional Conf. on Thermoelectrics, AIP Conf. Proc., 
316, 328 (1995). 

D.L. Partin, J. Electronic Materials, 13, 493 (1984). 
G. Dresselhaus, Proc. of ARPA Thermoelectric 
Materials Workshop, Arlington, VA, November 1995. 
R. People, J.C. Bean, and D.V. Lang, J. Vac. Sci. 
Technol., A3, 846 (1985). 

T. Harman, M.M. Manfra, and D.L. Spears, Submit- 
ted to J. Electron. Materials. 

R. Venkatasubramanian, T.S. Colpitts, M. Mantini 
and D. Malta, Proc. of 14th International Conf. on 
Thermoelectrics, St. Petersburg, in press (1996). 
U.S. Patent Provisional Application Serial 


No. 60/016, 710. 

A. Mzerd, D. Sayah, J.C. Tedenac, and A. Boyer, J. 
Cryst. Growth, 140, 365 (1994). 

R. Venkatasubramanian and T.S. Colpitts, Presented 
in the VIII International Conf. on OMVPE, Cardiff, 
England, June 1996. 

O. Lang, R. Schlaf, Y. Tomm, C. Pettenkofer, and W. 
Jaegermann, J. Appl. Phys., 75, 7805 (1994). 

B.E. Warren, X-Ray Diffraction , Reading, Mass., 
Addison-Wesley, (1969). 

J. Prater, Proc. of ARPA Thermoelectric Materials 
Workshop, Arlington, VA, November 1995. 

H. Scherrer and S. Scherrer, Proc. of the 12th 
International Conf. on Thermoelectrics, Yokohama, 
H. Mastubara Ed., 90 (1994). 

D.G. Cahill, M. Katiyar, and J.R. Abelson, Phys. 
Rev. B, 50, 6077 (1994). 

W.S. Capinski and H.J. Marris, “Thermal Conductiv- 
ity in GaAs/AIAs Superlattices,” Paper presented at 
the March 1996 American Physical Society Symp. on 
Semiconductors: Transport in Superlattices I; R. 
Venkatasubramanian, “Thermal Conductivity 
Reduction in Superlattice Structures,” Paper pre- 
sented at the March 1996 American Physical Society 
Symp. on Semiconductors: Transport in Superlattices 
I 

G. Chen, C.L. Tien, X. Wu, and J.S. Smith, J. Heat 
Transfer, 116, 325 (1994). 

Another approach is to process the 10-mm-thick 
films into thermoelements of proportionately smaller 
areas compared to bulk materials. 


Four/1996 = 41 





Profiles 
in Science 


Jean-Pierre Fleurial 


Dr. Fleurial earned his Ph.D. in solid state physics and 
materials science from the National Polytechnique Institute 
of Lorraine, France, in 1988, during which he developed new 
growth techniques and theoretical models for high quality 
single crystals of Bi,Te,. In 1988, Dr. Fleurial came to JPL 
as a National Research Council Resident Research Associ- 
ate and joined the permanent staff in 1990. He was first 
involved with the optimization of Si-Ge alloys for thermo- 
electric power generation. He has calculated ternary and 
quaternary phase diagrams of Si-Ge alloys, and also per- 
formed diffusion and solid solubility studies, liquid phase 
epitaxial growth of Si-Ge. More recently, Dr. Fleurial has 
been actively engaged in the study of novel promising ther- 
moelectric materials, in particular compounds with the 
skutterudite structure. He is also involved in thermoelectric 
microcoolers and high thermal conductivity materials for 
resolving thermal management issues. Dr. Fleurial currently 
leads the JPL Thermoelectric Team and is the active Secre- 
tary of the International Thermoelectric Society. 
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Thierry Caillat 


Dr. Caillat earned his Ph.D. in materials science from 
the National Polytechnique Institute of Lorraine, France, in 
1991. During his Ph.D., he conducted experimental and theo- 
retical research on high quality single crystals of p-type 
Bi,Te,-Sb,Te, alloys grown by a novel technique, the travel- 
ing heater method. After this research, Dr. Caillat received a 
National Research Council Fellowship to pursue studies of 
new semiconducting materials for thermoelectric applications 
at JPL. He joined the permanent staff at JPL in 1994. In the 
last five years, he has played a key role in identifying several 
families of promising compounds for thermoelectric appli- 
cations, including skutterudites and b-Zn,Sb,-based materi- 
als. He has been elected to serve as a Board Member of the 
International Thermoelectric Society from 1996 to 1999. 





Profiles in Science 
(continued) 


JPL’s Thermoelectric Team, Materials and Device 


Technologies Group 


The staff of Jet Propulsion Lab's (JPL’s) 
Thermoelectrics Team in the Electric Power Systems Sec- 
tion is considered by the international community as a leader 
in developing thermoelectric materials and understanding 
thermoelectric power generation and cooling technology. 
Jean-Pierre Fleurial and Thierry Caillat are noted members 
of the team receiving support from the Office of Naval Re- 
search (ONR). The Team has a proven record of establishing 
unique facilities as dictated by the needs of a particular study 
and is widely known for its outstanding materials prepara- 
tion and characterization capabilities. A wide variety of prepa- 
ration techniques are exploited: crystal growth and powder 
metallurgy for bulk materials and various deposition pro- 
cesses for thin films. The electrical resistivity, Hall effect, 
Seebeck coefficient, thermal diffusivity and the heat capac- 
ity can be measured from 4K to 1300K. In addition, equip- 
ment built for testing devices at high temperatures and for 
extended periods of time is considered a key capability. 

The Team has also pioneered the development of a 
comprehensive solid state physics model for the thermal and 
electrical transport properties of a given material over its full 
temperature range of usefulness to use as a tool for guiding 
experimental and for predicting the maximum performance 
likely to be achieved. In the past 5 years, several programs 
were initiated on the development of new thermoelectric 
materials with high thermoelectric efficiencies for space and 
terrestrial applications. Recent results obtained at JPL and 
confirmed by other laboratories have shown that ZT (figure 
of merit) values significantly larger than | could be obtained. 


A new growing area for the Team is the development of ther- 
moelectric microcoolers and high thermal conductivity ma- 
terials such as diamond for the thermal management of high 
density power electronics. These JPL programs receive sup- 
port from ONR, National Aeronautics and Space Adminis- 
tration (NASA), other federal agencies, and private indus- 
tries. 
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Thermoelectric 


Property 


Measurements 


Ctirad Uher, Department of Physics, University of Michigan, Ann Arbor, MI 


As discussed elsewhere in this journal, the most im- 
portant parameter determining whether a given solid is a po- 
tentially useful thermoelectric material is the figure of merit 
Z = S’o/« or its dimensional equivalent ZT. Here S is the 
thermoelectric power or the Seebeck coefficient, o is the elec- 
trical conductivity, « is the thermal conductivity and T is the 
absolute temperature. Thus, to assess a given material for its 
thermoelectric potential necessitates measuring the three 
essential parameters S, 6, and x. The measurements should 
be carried out over a wide enough temperature range to gain 
insight about the trend in the transport behavior and to es- 
tablish the domain over which the figure of merit is maxi- 
mized. Often one desires to have more detailed information 
about the transport properties and an inquiry into the carrier 
density, the dominant carrier type, mobilities, effective masses 
and scattering times is relevant and helpful in revealing the 
overall transport behavior. Such expanded property measure- 
ments then include the whole gamut of transport investiga- 
tions with the Hall effect as clearly the most important one. 
The singular importance of the Hall effect is the reason why 
I include it in the review of the transport properties which 
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determine the figure of merit of the thermvelectric materi- 
als. 


l. 


Electrical Kesistivity 


Electrical resistivity is perhaps the most frequently 
performed measurement on a conducting solid. Because of 
its conceptual simplicity such a measurement might appear 
as a straightforward exercise requiring no special consider- 
ations. This is not so in practice and, especially in investiga- 
tions of thermoelectric materials, great care must be used to 
assure accurate and reliable experimental results. 

The measurements of the resistivity can be done by a 
variety of techniques and the preferred one depends on fac- 
tors such as the size of a specimen, its shape, magnitude of 
the resistivity, desired accuracy of the data, available instru- 
mentation and, to a degree, taste and preference of the per- 
son conducting the measurements. Before I describe the ex- 
perimental techniques most frequently applied to measure- 
ments of the resistivity it is worthwhile to note the range of 





resistivities one encounters in thermoelectric materials. As 
is stated elsewhere in this journal, the carrier concentration 
in practical thermoelectrics falls within the range 10*°-107° 
m* which implies that these materials have rather low val- 
ues of resistivity somewhere in the 10° Q-m range. Thus, 
accurate measurements necessitate the use of high sensitiv- 
ity instrumentation. Furthermore, because the resistivity is 
low, the contact resistance is comparable to the sample re- 
sistance and one should avoid using a two probe configura- 
tion where the same two contacts serve the dual role of car- 
rying the sample current and sensing the potential difference. 
For any kind of reliable measurements the use of four probe 
techniques with well separated current and voltage contacts 
is essential. 


la. Four-probe bar 
measurements 


Perhaps the most frequently used technique employs a 
sample in the form of a parallelepiped (bar) with two pairs 
of contacts (four probes), Fig.la. One pair serves to inject 
current through large-area ohmic contacts at the ends of the 
bar and the other pair, separated from the current contacts, 
measures the voltage drop along the length of the bar via 
small-area contacts. It can be shown by a relatively straight- 
forward analysis’ that by insetting the voltage contacts by a 
distance not shorter than the width of the sample, the sample 
current over the distance delineated by the voltage probes 
will be axial and the equipotentials will remain essentially 
planar and normal to the bar axis even for highly nonuni- 
form current contacts. The voltage probes are typically held 
in place by either pressure or by means of a solder or a vari- 
ety of conducting pastes and epoxies. 

Frequently, point contacts may be too resistive and it 
is preferable to prepare the specimen in the form of a bridge 
with long and narrow side-arms to which large-area contacts 





Figure 1a. 
Bar-shaped sample with a typical arrangement of probes for 
measurements of the electrical resistivity. 


Figure 1b. 

Bridge-shaped specimen suitable for measurement of the 
resistivity and Hall effect. Pads for attachment of contacts 
are shaded. 








can conveniently be attached. An example of a bridge-shaped 
sample is given in Fig. 1b. Although no current flows through 
the side-arms they do cause distortions of the equipotentials 
in their vicinity. Fortunately the effect is rather small for side- 
arms with the width substantially smaller than the width of 
the sample and one should aim to make them as narrow as 
mechanical strength allows. Bridge-shaped samples are a 
popular geometry and they have an added advantage of be- 
ing convenient for measuring of the Hall effect in which case 
either one of the opposing pairs of side-arms could be used 
to pick up the transverse (Hall) voltage. 

The electrical resistivity p for the above two sample 
geometries is given as 


A 
oe ee 


ial a aa 


(1) 

where R is the resistance, L the distance between the voltage 
probes and A the cross-sectional area of the sample. Thus, 
this method requires an accurate electrical measurement to 
determine R as well as a precision mechanical measurement 
to specify the dimensions of the sample (geometrical factor 
A/L). A typical experimental set-up to measure the resistiv- 
ity of a bar or bridge-shaped sample includes a standard re- 
sistor of resistance R, in series with the sample in order to 
determine the actual current in the circuit 1 = V/R,. The 
simplicity of this experimental technique is its main draw 
and appeal. However, when used to measure the electrical 
resistivity of thermoelectric materials, one needs to be aware 
of its shortcomings. The chief problem is that the voltage 
probes pick up the total potential which may or may not con- 
sist of contributions from sources other than Ohm’s law and 
it is difficult to account for some of these effects. If one could 
assume strictly isothermal conditions there would be no prob- 
lem. Unfortunately one frequently encounters standing ther- 
mal gradients AT across the sample (either imposed by im- 
proper experimental set-up or due to excessive Joule heat- 
ing) and, because of the large thermopower of thermoelec- 
tric materials, these standing thermal gradients contribute 
thermoelectric voltage AV = S A T which is superposed on 
the truly resistive voltage drop. One can deal with the ther- 
moelectric contributions which originate in the standing ther- 
mal gradients simply by reversing the current flow and tak- 
ing the average of the two readings. A far more serious prob- 
lem concerns the Peltier effect contribution. As is well known, 
a passage of current through a junction between two dis- 
similar metals is accompanied by the heat flux density and 
this leads, depending on the direction of the current, to lib- 
eration or absorption of heat at the junction. Since one usu- 
ally does not have available fine wire-probes made of the 
same material as the sample under study, one obviously faces 
the problem associated with the Peltier heat at the two cur- 
rent contacts. Because the Peltier effect is linear in current, 
current reversal changes the sign of both the resistive volt- 
age and the Peltier effect and, one cannot separate the re- 
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spective contributions by averaging the two readings. The 
net effect is thus an apparent increase in resistivity which, 
for a good thermoelectric material such as Bi,Te, amounts 
up to 1/3 of the total resistivity. 

The most effective way to deal with the Peltier contri- 
bution is to use alternating current in combination with a 
lock-in amplifier or to introduce an electromechanical chop- 
per which converts the de current into a square-wave alter- 
nating current of suitable frequency. Since it takes a finite 
time to establish the thermal gradient due to the Peltier heat 
while the resistive voltage sets in instantaneously, the idea is 
to switch the current polarity before the thermoelectric con- 
tribution can build up. Because one can use accurate poten- 
tiometers or high sensitivity digital voltmeters to pick up the 
demodulated signal and one avoids the problems of stray 
signals associated with wide frequency band ac meters, the 
chopped dc is preferable to alternating current. Chopping 
frequencies of several tens of Hz are typically perfectly ad- 
equate to suppress the Peltier effect contribution. Potentio- 
metric circuits utilizing the electromechanical choppers have 
been described in the literature’. 


Ib. The van der Pauw 
technique 


The requirement of a bar-shaped or bridge-shaped 
sample to ensure uniform current flow and well defined 
equipotentials represents a serious constraint on the sample 
geometry. It is often very inconvenient to process a sample 
into a bar-shaped form and one would rather use a measur- 
ing technique which is less restrictive in terms of the required 
sample geometry and suitable for samples of irregular con- 
tour or plane lamella. The technique which satisfies these 
criteria is based on the conformal mapping of two-dimen- 
sional fields and was developed by van der Pauw’. This tech- 
nique has found widespread use in recent years and its only 
essential requirement is that the sample be of uniform thick- 
ness, structurally homogeneous, and singly connected, i.e., 
not having isolated holes. 

Consider a plane lamella of uniform resistivity p and 
thickness 5 with four line contacts on the edge separated by 
distances a, b, and c as shown in the inset of Fig. 2. We are 
interested in the potential at terminals 3 and 4 as a result of 
the current I introduced at the terminal | and removed at the 


terminal 2. It is easy to show that the resistance R,, ,, defined 


a 
as the ratio —* equals 
12 


V, 


b)\(b +c) 
R =—t= Py (a+ b\b+c) 
a ae . 


2 
b(a+b+c) “ 


By commuting the electrodes, i.e., using terminals | and 4 
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as voltage electrodes while injecting the current into the half- 
plane via the terminal 2 and removing it at the terminal 3, 
the resistance R,, ,, becomes 


Pp», (a+b)(b+c) 
md ac 


(3) 


Ry = 


Ym 
1, 


A simple manipulation and the fact that b (a+b+c) + ac = 
(a+b) (b+c) allows one to combine Egs.2 and 3 into a single 
equation containing the desired resistivity p of the lamella, 


ze A + Ry 41 Kt Rip 34 
P= in2 2 
n Ry3.41 


Here f is the asymmetry function which depends only on the 
ratio Ris/RasarA plot of the asymmetry function f versus 
the resistance ratio Ri s/Ras41 is given in Fig.2. The experi- 
mental procedure to determine the resistivity using the van 
der Pauw technique thus consists of measuring the two 
resistivites R,, ,, and R,, ,,, using Fig.2 to determine the value 
of the asymmetry function f corresponding to the resistance 
ratio Ri s/Ra,, a1, and inserting the values into Eq.4 to deter- 
mine the resistivity of the material. 





(4) 





Figure 2. 


The asymmetry function f as a function of the resistance 
ratio R,, ,/R,, 4, The inset shows a sample in the form of a 
semi-infinite slab of thickness 6 with four line-contacts 
separated by distances a,b and c along the edge of the slab. 
Current enters the sample at the contact 1 and leaves it at 
the contact 2. Potential difference is measured between 
contacts 3 and 4. 
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While Eq.4 applies to an arbitrarily shaped lamella, 
there is an advantage in using symmetrical shapes such as 
discs or squares contacted in a symmetrical fashion. For in- 
stance, if the specimen is a perfect circle with the contacts 
along the orthogonal diameters, then obviously R,, ,,=R,, ,, 
the asymmetry function f has value of unity, and Eq 4 be- 
comes 


8 


R = In2 Riss (5) 


where R stands for the sheet resistance measured in “ohms 
per square.” 

The location and size of the contacts is crucial for the 
accuracy of resistivity measured by the van der Pauw tech- 
nique. A theoretical assumption of line-contacts attached at 
the edges of the lamella is not always viable in practice and 
it is preferable to place the contacts on top of the structure 
close to the outer edge. The price one pays for violation of 
the assumed edge-placement of the contacts and for their 
finite size is an error introduced in the measurement of resis- 
tivity. Van der Pauw and subsequently a number of authors*° 
considered various electrode layouts and sample shapes and 
calculated the error introduced by each specific configura- 
tion, including the so-called clover-leaf sample geometry. 
In general, the finite size of contacts and their displacement 
from the edge tend to result in the measured resistivity which 
is smaller than the actual resistivity of a material under study. 
Planar shape of the samples including thin-films deposited 
on a variety of substrates tend to smooth out any standing 
thermal gradients better than do the bar-shaped specimens. 
The van der Pauw technique is thus a useful method of as- 
sessing the electrical resistivity of thermoelectric materials. 

One should note that the electrical resistivity is also 
often measured by the collinear four-point probe technique”*. 
Although used rather sporadically on thermoelectric materi- 
als (due to small voltages developed between the closely 
spaced voltage probes as a result of relatively low resistivity 
of thermoelectric materials), numerous versions of this tech- 
nique, including a number of commercially available set- 
ups, have found widespread use in evaluating materials which 
typically have resistivities in the range 10° - 10' Q-m. Al- 
though not as accurate as the techniques where the geometri- 
cal factor of the sample is well defined, the four-point probe 
might be useful for rapid scanning of production ingots of 
thermoelectric materials. For a detailed description of this 
deceptively simple technique and how to account for vari- 
ous correction factors the reader is referred to the original 
literature*"'. 


2. Seebeck Coefficient 


As the definition S = AV/AT implies, measurement of 
the Seebeck coefficient, or the thermoelectric power, as the 


Seebeck coefficient is frequently called, involves a simulta- 
neous determination of a thermoelectric voltage AV arising 
as a consequence of an imposed temperature difference AT 
between two points along the sample length. To get an accu- 
rate value of S, both AV and AT must be measured with pre- 
cision and that requires certain special considerations. 
Experimental arrangements to determine the Seebeck 
coefficient usually involve one of two configurations illus- 
trated in Fig. 3. In the first one, a sample is provided with a 
heater (a small resistor, wire-wound heater, or thin-film strain 
gauge resistor) at one end and the other end is attached (sol- 
dered, epoxied or by other means that provide good thermal 
link) to a heat sink. Thermal and potential differences are 
then measured between two points along the length of the 
sample. Depending on the temperature range of interest, 
sample size and desired accuracy, one can use small cali- 
brated resistance thermometers or, more typically, thermo- 
couples to determine AT, and either fine gauge copper wires 
or the legs of the thermocouples to measure the Seebeck 
voltage. In the other popular configuration, shown in Fig. 
3b, the sample is clamped between two reasonably large cop- 
per blocks and both AV and AT are measured between the 
blocks. Because of the inherently low thermal conductivity 
and high Seebeck coefficient of the thermoelectric materials 
and consequently an insignificant contribution of interfaces 
to the overall thermal and potential gradients, it makes little 
difference which of the two experimental set-ups is employed. 
For other materials where the contact resistance might intro- 
duce large interfacial gradients the layout of Fig. 3a is al- 
ways the preferred one. It should be noted that measurements 
of the Seebeck coefficient do not require good vacuum. In 
fact, the presence of some inert gas might be beneficial to 





Figure 3. 


Experimental set-up for the measurement of the Seebeck 
coefficient using copper-constantan thermocouples. 

3a. Temperature sensors attached to the sample. 

3b. Temperature sensors attached to the heat source and 
sink. 

Numbers 1 and 2 designate copper and constantan wires, 
respectively. 
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improving the thermal contact. However, because it is con- 
venient to carry out measurements of the Seebeck coeffi- 
cient together with the measurements of the thermal con- 
ductivity, and the latter requires good vacuum, the studies 
are most often done under adiabatic conditions. 

Assuming one uses temperature sensors of reputable 
origin, the success of the experiment hinges primarily on the 
quality of the thermal contact at the thermocouple junctions. 
An ideal thermal contact is achieved by soldering the ther- 
mocouple junction to the sample. If one desires to use a dif- 
ferential thermocouple, it is essential to remember that at 
least one of the two junctions must be insulated from the 
sample otherwise the emf generated by the thermocouple is 
shorted via the sample and an error reading results. Alterna- 
tive ways of attaching the thermocouple junctions rely on 
various glues and epoxies and one may also opt for inserting 
the thermocouple junctions into tiny holes drilled in a speci- 
men. 

Attention paid to the reference junction is just as im- 
portant. For room temperature measurements and for inves- 
tigations at elevated temperatures one usually keeps the ref- 
erence junctions at an ambient temperature. All that is needed 
is a small Dewar flask filled with water and covered with a 
polystyrene cap through which is inserted a thermometer and 
glass tubes filled with oil. Thermocouple junctions are im- 
mersed in the oil which facilitates good and stable thermal 
contact. Alternatively, one can use commercially available 
ice reference points which are capable of maintaining the 
temperature at 0+ 0.1°C. Measurements below the room tem- 
perature are normally done in a cryostat an integral part of 
which is a copper sample block or copper “pot” which can 
serve as a robust thermal grounding post. The temperature 
of this copper post is carefully monitored using a platinum 
thermometer or some other high sensitivity resistance ther- 
mometer. Because of their low sensitivity at liquid helium 
temperatures, thermocouples must eventually be replaced by 
resistance thermometers if one wishes to study the Seebeck 
effect below about 10K. This of course avoids the problem 
with thermocouple junctions but adds a new one - how to 
cope with much more bulky resistance thermometers espe- 
cially if the sample is small. 

The choice of a thermocouple depends very much on 
the temperature range of interest and on how much stray heat 
can one tolerate to reach the sample. Basically, the desire is 
to maximize the thermocouple reading and to minimize the 
heat path into/from the sample along the leads. Because of 
its very high thermal conductivity, if copper forms one leg 
of a thermocouple, it should be a very thin wire. In general, 
one should only use thermocouple wires that come from a 
reputable source and for which exist calibration tables trace- 
able to the national standard.'* One should avoid using old, 
strained and “kinking” thermocouples as these may gener- 
ate spurious signals. If one makes his/her own thermocouple 
junctions, one should always test the efficacy of the couple 
by comparing the emf generated at two fixed temperatures 
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(e.g., the ice point and liquid nitrogen) and compare the read- 
ings with the calibration tables. Special attention needs to be 
paid to measurements at high temperatures. One should not 
only periodically recheck the thermocouple calibration but 
it is also essential that the wires brought out from a high 
temperature environment to the measuring instrumentation 
are not interrupted at points where large temperature differ- 
ences are likely to exist. It is always preferable to use capil- 
lary tubes and pass the wires uninterrupted. 

Although the propensity for errors arising from mea- 
surements of AT is far greater than for measurements of AV, 
and the errors are usually more serious , it is, nevertheless, a 
prudent practice to ascertain that AV is the true Seebeck volt- 
age. The key point to keep in mind is the fact that the Seebeck 
voltage measured is the property of a couple formed by the 
sample and the leads, Sous = Scope ~ Sieass” 4S We already 
noted, the connecting leads ‘could be the legs of the thermo- 
couples or a pair of separate wires. In either case, one must 
correct for the absolute thermopower of the leads in order to 
obtain the absolute thermopower of the sample. 
Thermopower of the copper wires, the material most fre- 
quently used as the Seebeck voltage probes, does not exceed 
a couple of p'V/K at 300K while typical values of the Seebeck 
coefficient of thermoelectric materials are in excess of 200 
utV/K. Thus, for thermoelectrics, the omission to correct for 
the thermopower of copper leads does not lead to a serious 
error. This is not so for other materials and in any case it is a 
good practice to present the data in their correct form. Tables 
providing thermopower values as a function of temperature 
for pure Pb wire which serves as the standard reference ma- 
terial for calibration purposes are available in the literature.’ 
Furthermore, the discovery of high temperature supercon- 
ducting perovskites (superconductors have zero absolute 
Seebeck coefficient at temperatures T< T_) has greatly ex- 
tended the temperature range over which the absolute 
Seebeck coefficient of any material can quickly and conve- 
niently be determined." 

I wish to note that the Seebeck coefficient ean also be 
determined using the Harman method discussed in section 
4. In this case a passage of dc current through a sample sus- 
pended in a vacuum (adiabatic conditions) results in the volt- 
age across the sample V_, and the temperature difference AT. 
By subtracting the ohmic voltage V, (measured under the 
isothermal conditions) from the voltage V, and dividing by 
AT one obtains the Seebeck coefficient of the sample with 
respect to the measuring leads. Correcting for the 
thermopower of the leads yields the absolute Seebeck coef- 
ficient of the sample. Though simple in principle, it is not 
advisable to rely on this method for an accurate determina- 
tion of the Seebeck coefficient. The reason is a small tem- 
perature difference AT and the fact that the adiabatic and 
isothermal voltages V, and V, do not always (especially at 
lower temperatures) differ by a margin large enough to 
achieve, after subtraction, good enough accuracy. 





3. 


Thermal Conductivity 


Of the three transport properties that determine the fig- 
ure of merit of a thermoelectric material, the thermal con- 
ductivity is technically the most demanding one but also the 
most sensitive to the nuances of microstructure and compo- 
sition. This latter attribute has in fact been used to great ef- 
fect in increasing the figure of merit of a given material sys- 
tem. Traditionally, thermoelectric materials have been de- 
veloped by selecting systems with reasonably large power 
factors (the power factor is the numerator So in the figure 
of merit) and then systematically decreasing the thermal con- 
ductivity by altering the microstructure, defect and impurity 
content, or chemical composition. This means that when 
searching for potentially useful thermoelectric materials or 
when optimizing their thermal conductivity, the experimen- 
tal technique must be flexible enough to accommodate 
samples with a wide range of thermal conductivities. While 
the heat flux, thermal gradient, and thermal conductivity are 
analogous to electrical current flow, potential difference, and 
electrical conductivity in an Ohm’s Law experiment, accu- 
rate thermal measurements are far more difficult to realize 
in practice than their electrical analogue. This is chiefly due 
to the fact that while in an electrical resistivity experiment 
current leakage paths can be virtually eliminated, it is im- 
possible to completely suppress all heat leakage paths in a 
thermal conductivity experiment. The diligent experimentalist 
must always be aware of the existence of heat leaks in his 
experimental apparatus and take special care to minimize 
them with respect to the heat flow through the sample. 

I now discuss several methods for measuring the ther- 
mal conductivity of thermoelectric materials and highlight 
their relative advantages and disadvantages. The choice of 
the appropriate method is a strong function of the sample 
size and geometry, material type, and temperature range of 
interest, and in many cases it may be desirable to use more 
than one method to insure an accurate determination of the 
absolute magnitude of the thermal conductivity. 


3a) The Longitudinal Steady- 
State (LSS) Technique 


This technique is the thermal analogue of the four- 
probe bar measurement of electrical resistivity discussed in 
section la. For this type of measurement, a long bar-shaped 
or cylindrical sample is preferred. It has the advantages of 
providing an absolute measurement of the thermal conduc- 
tivity over a wide temperature range using fairly simple and 
inexpensive instrumentation. Its chief disadvantages are limi- 
tations on sample size (generally only samples longer than a 
few millimeters can be measured) and sample geometry (the 
method is unsuitable for thin film materials and for samples 
with highly irregular shape), and the susceptibility to errors 


due to heat losses especially in measurements on low ther- 
mal conductivity materials. The generic experimental set-up 
is shown in Fig.4. One end of the sample, of cross-sectional 
area A, is secured by means of epoxy, solder, or by clamping 
to a temperature controlled heat sink. This heat sink can be 
the cold tip of a liquid nitrogen or liquid helium cryostat, the 
cold head of a closed cycle refrigerator, or even a heated 
plate. An electrical heater is attached in a similar manner to 
the free end of the sample. This heater may be a small wire 
wound resistor, a strain gauge, or an evaporated metal film. 
For a given current I and voltage V across this electrical 
heater, an amount of power P = IV will be produced. At two 
points, spatially separated by a distance L along the sample, 
are affixed temperature sensors which measure the tempera- 
tures T, and T_ of the hot and cold side of the sample. When 
the heater is energized, a heat flux Q=P/A will flow down 
the sample and a temperature difference AT will be devel- 
oped. Under steady-state condition, and assuming all of the 
Joule heat generated by the heater flows down the sample, 
the thermal conductivity can be written as 


(6) 


Thus by measuring the heater power, the temperature differ- 
ence along the sample, and geometrical factor L/A the ther- 
mal conductivity is determined. Clearly the accuracy of this 
method depends on the accuracy to which the quantities P, 
AT, and L/A can be measured. The power P is specified by 





Figure 4. 


Longitudinal steady-state technique to measure thermal 
conductivity of a bar shaped sample of cross-sectional area 
A. After the heater is energized with power P, a temperature 
difference AT = T, - T, is measured between two points 
separated by a distance L along the sample. 
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the product of I and V and both quantities can be measured 
to within one part in 10*. The uncertainty in the geometrical 
factor L/A can be as small as | % for large, long samples 
and as large as 20 % for short samples or specimens with 
irregular shape. 

The most critical parameter in the steady state mea- 
surement of « is the temperature difference AT. The accu- 
racy to which this quantity can be measured is a function of 
the type of temperature sensors used and the temperature 
range of the measurement. The two most popular types of 
sensors for measuring sample temperatures in a steady-state 
experiment are resistance thermometers and thermocouples. 
Each of these sensors has its advantages and drawbacks and 
the choice largely depends on sample size and geometry and 
temperature range of interest. 

Resistance thermometers offer a highly accurate mea- 
surement of the absolute temperature because they possess a 
strongly temperature dependent and well characterized and 
reproducible resistance. They are commercially available 
from many sources'* and can be purchased already calibrated, 
so that a measurement of the resistance of the sensor pro- 
vides a very accurate determination of the absolute tempera- 
ture. A notable drawback of resistance thermometers is their 
relatively bulky size which makes it difficult to accommo- 
date these thermometers on samples of short length. 

For smaller samples the use of fine wire thermocouples 
is the preferred technique and the issues pertaining to their 
proper use as discussed in section 2 are equally relevant here. 
Caution should also be exercised when using thermocouples 
in a magnetic field. The Seebeck coefficients of the thermo- 
couple materials can be a quite strong function of field and 
their field dependence can also differ from strand to strand. 
The experimentalist who wishes to measure the thermal con- 
ductivity of a specimen in a magnetic field using thermo- 
couples should take every precaution to ensure the field de- 
pendence of the sensor signal is well characterized. 

Up to now I have tacitly assumed that all of the heat 
generated by the heater flows down the length of the sample 
and into the heat sink; in an actual experiment this condition 
is never achieved. The challenge to the experimenter is to 
minimize the alternate means of heat transfer to ensure that 
virtually all the heat flows down the length of the specimen. 
Let us consider the main sources of heat loss. 

When the heater in Fig.4 is energized, heat has the 
following alternate paths: i) conduction down the length of 
the sample; ii) conduction and convection in the surround- 
ing medium; iii) conduction along the heater and tempera- 
ture sensor leads; and iv) radiation to the surroundings. The 
goal is to minimize mechanisms ii)-iv) with respect to the 
desired mechanism i). Heat transfer via mechanism ii) can 
be virtually eliminated by placing the measuring apparatus 
in a vacuum chamber evacuated to less than 10° Torr. Con- 
duction along the sensor and heater leads can be made insig- 
nificantly small by using long, thin wires which can be coiled 
to fit into the experimental space. The chief limitation on the 
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steady-state measurement is heat loss from the heater and 
the sample by radiation. These radiative losses depend on 
the 4th power of the absolute temperature, and thus at low 
temperatures (say below 100 K or so) they can generally be 
neglected. At higher temperatures and for samples with small 
thermal conductance, radiation losses can become substan- 
tial. These losses can to some extent be minimized by sur- 
rounding the sample with a radiation shield thermally an- 
chored to the heat sink, or better yet, provided with a tem- 
perature profile along its length matching that across the 
sample. This will minimize the temperature difference be- 
tween the radiating surfaces and thereby minimize heat trans- 
fer by radiation. In any case it is desirable to have some esti- 
mate of the magnitude of the radiation losses in a given ex- 
perimertal setup. This can be done either by modeling the 
system using Stefan’s Law and the appropriate value of emis- 
sivity, or by measuring the total heat losses directly, e.g., by 
replacing the sample with a very long, thin sample of low 
thermal conductivity such as glass, energizing the heater, and 
observing the temperature profile along the sample. Even 
with all these precautions, one really cannot hope to make 
accurate measurements of thermal conductivity using the 
steady-state technique much above room temperature on any 
material of thermal conductivity less than a few hundred mW 
cm K". Thus in order to characterize materials for high 
temperature applications, other techniques are called for. 


3b) The Radial Steady-State 
(KSS) Method 


Among the best known experiments utilizing the ra- 
dial steady-state method are those measuring the thermal 
conductivity of silicon and germanium up to the melting 
point.'* In this technique, the sample is a cylinder of length 
L and heat is made to flow radially out from the center of the 
sample; see Fig.5. Thus a heater wire must be placed along 
the axis of the sample and two thermocouples placed at two 
radial positions r, and r,. In this geometry the thermal con- 
ductivity is given by 


_ Pin(y,/7) 
~ QnLAT 


(7) 


The sample must be slit in half in order to insert the heater 
and thermocouple wires and then clamped tightly back to- 
gether. The sawn faces must be smoothly polished so that 
they fit neatly together in order to eliminate any conduction, 
convection, or radiation along the interface. Although this 
technique is conceptually very simple and elegant, it is diffi- 
cult to realize in practice because large samples (for instance 
1 cm in diameter by several centimeters long) are required. 





3c) The 3-® Method 


While the LSS technique described in section 3a of- 
fers the distinct advantages of relative simplicity and a di- 
rect measurement of the thermal conductivity, it has two main 
drawbacks: radiation losses limit its use at high tempera- 
tures, and in general it cannot be used to measure thin film 
materials. A technique developed over the last ten years which 
can overcome both of these difficulties is the 3-@ technique.'” 
Intended originally for measurements of the thermal con- 
ductivity of bulk glasses in the 30-750 K temperature range, 
it has since been refined and adapted to measure thermal 
conductivity of thin films. Typical samples are 1 cm square 
plates a few millimeters thick, or a thin film deposited on a 
substrate. Onto the surface of the sample is evaporated a thin 
metal line (such as silver or platinum) which serves the dual 
purpose of heater and thermometer. A current through the 
heater at frequency @ produces a temperature oscillation at a 
frequency of 2a. Since the resistance of the heater-thermom- 
eter increases with temperature it also oscillates at 2. This 
resistance oscillation multiplied by the original current os- 
cillation at frequency @ produces a voltage oscillation across 
the strip at frequency 3@, whence the method’s name. By 
solving the heat transfer equation for a narrow line source 
on the surface of an infinite half volume, this voltage can be 
directly related to the thermal conductivity of the material. 
A nice feature of the method is that by increasing @, regions 
closer and closer to the surface of the sample can be probed, 





Figure 5. 


A radial steady-state thermal conductivity measurement. A 
heater H is placed along the axis of a cylindrical sample of 
length L. A temperature difference AT = T, - T, is measured 
between two points r, and r, along the radius of the cylinder. 








i.e., at high enough frequency the method will measure the 
thermai conductivity of a film on a substrate. Before apply- 
ing this technique to a broad range of materials, an impor- 
tant point must be remembered: the technique as originally 
developed was to measure thermal conductivity of electri- 
cally insulating materials. If one wishes to apply this tech- 
nique to a thermoelectric material (a rather good electrical 
conductor), an electrically insulating layer must be placed 
between the heater/thermometer and the sample. The experi- 
ment then measures the sum of the thermal resistance of the 
insulating layer and the sample, as well as the boundary re- 
sistances at each interface. In principle it should be possible 
to separate out these contributions by varying the frequency 
or simply minimizing the contribution of the insulating layer 
by making it very thin, but such modifications have not yet 
been investigated in detail. The 3-@ technique does, how- 
ever, provide the potential of directly measuring the thermal 
conductivity of thin film thermoelectric materials. Films as 
thin as 10 microns can be measured using the original 3-@ 
technique without modification, and recently it has been dem- 
onstrated'* that the method can be extended to probe films as 
thin as 3000 A. In addition, the method is essentially free 
from errors due to radiation losses, and thus may be useful 
in measuring bulk materials at high temperatures. The rea- 
son the radiation losses do not enter into consideration is 
because they are proportional to the characteristic length scale 
of the measurement, which can be made very small in the 3- 
@ technique by using photolithography to define the heater/ 
thermometer. 


3d) Transient Method of 
Measurement 


Several transient techniques have been described in the 
literature'®”° and all are variations to some degree of the origi- 
nal transient thermal wave experiments of Angstrom.?! These 
methods monitor the temperature rise on the back surface of 
a sample when the front surface is exposed to a heat wave 
that has a time dependence. A modern, laser-based version 
of the method” consists of a pulsed laser heat source, a 
vacuum furnace, and detection, amplification and analysis 
systems. The front and back surfaces of the sample must be 
covered with a coating to absorb the radiant incident energy. 
The laser is pulsed with a duration on the order of 1 ms and 
total energy output on the order of 10 J. The temporal de- 
pendence of the temperature rise on the back surface of the 
sample as recorded by a suitable infrared detector is related 
to the thermal diffusivity of the sample. As in the case of the 
3-@ technique, radiation losses can be kept small by reduc- 
ing the size and thickness of the sample; this however, re- 
quires a pulse of shorter duration since the pulse duration 
must remain small with respect to the transit time of the heat 
wave through the sample. A potential pitfall of this technique 
is the reliable and reproducible production of an absorbing 
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layer on the sample, which can cause difficulty when adhe- 
sion and thermal expansion problems arise at high tempera- 
tures. In addition, the determination of the thermal conduc- 
tivity from the measured diffusivity requires knowledge of 
the heat capacity C of the material since k = DCp, where p 
is the mass density. Thus C must be measured independently 
or calculated from the crystal structure and elastic constants. 
In spite of these drawbacks and the rather more expensive 
instrumentation required compared to traditional steady-state 
methods, transient thermal diffusivity measurements have 
become a popular approach of determining the thermal con- 
ductivity of thermoelectric materials above room tempera- 
ture. 

We conclude this section by describing two experimen- 
tal approaches used to identify the mechanisms by which 
heat flows in a thermoelectric material, information which 
is crucial when developing new thermoelectrics and opti- 
mizing their properties. 


3e) Separating the Electronic 
and Lattice Thermal 
Conductivities 


The thermal conductivity of a solid is composed of 
two contributions: that due to the charge carriers, K_, and 
that due to lattice vibrations, or, phonons, K.. The measured 
thermal conductivity is the sum kK = K_ + K.. For carrier con- 
centrations below about 10° m® the electronic component 
of thermal conductivity is usually negligible in comparison 
to K,. However, in some materials x, itself is very small (this 
of course is just what one desires for a thermoelectric), and 
in these cases kK, can be an appreciable fraction of the total 
thermal conductivity even at carrier concentrations below 
10° m®. In the course of understanding and optimizing the 
properties of a given material, it is essential to know the in- 
dividual contributions of phonons and carriers to the total 
heat conductivity. Generally this is done by either estima- 
tion of the carrier component using the Wiedemann-Franz 
(W-F) law or by suppression of k_ using a magnetic field. 

An estimation of the carrier thermal conductivity us- 
ing the W-F law is by far the easier of the two methods as the 
only required information is the magnitude of the electrical 
resistivity, and this is normally obtained as part of the rou- 
tine of evaluating the thermoelectric properties such as de- 
scribed in section 1. The W-F law states that, if the carriers 
in the system are undergoing scattering processes in which 
the change of energy of the carrier is small compared to the 
thermal energy «,T (elastic scattering), then x, = L, T/p 
where L, = 2.45 x 10* V’ K* is the Sommerfeld value of the 
Lorenz number and p the electrical resistivity. When the car- 
rier scattering is not elastic, the Lorenz number is smaller 
than L, and the Wiedemann-Franz law provides an upper 
limit to the carrier thermal conductivity. 
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In circumstances when the W-F law is known to break 
down and an exact determination of K_ is desired, one can 
attempt a separation of the thermal conductivity into its two 
components by using a magnetic field. A magnetic field ap- 
plied transverse to the thermal current can suppress the car- 
rier thermal conductivity when the field is large enough to 
achieve the classical high field regime specified by @.t >> 1 
or equivalently, JB >> 1. Here @. is the cyclotron frequency, 
t the scattering time and J the carrier mobility. Physically 
this condition implies that the carriers are “trapped” in sev- 
eral cyclotron orbits before being scattered and cannot trans- 
port their heat content down the thermal gradient. Conse- 
quently the heat conduction proceeds only via lattice vibra- 
tions. The high-field condition can actually be achieved only 
in fairly high mobility materials. For instance for B = 5 T, 
the product JB equals unity for p = 2000 cm? V's". Asa 
function of magnetic field one monitors a decrease and even- 
tual saturation of the thermal conductivity at a value equal to 
Ky Even if the mobility is not large enough to achieve the 
high field condition, extrapolation techniques have been de- 
veloped which allow x, to be determined.”** 


4. Direct Measurements 


of ZT 


Instead of measuring the transport parameters that 
specify the figure of merit in separate experiments, Harman™ 
showed that the figure of merit can be conveniently deter- 
mined in a single experimental run, the technique nowadays 
known as the Harman method. In section la we drew atten- 
tion to the existence of the Peltier heat that necessarily arises 
when a dc electric current passes through the sample and 
what complications this leads to in measurements of the re- 
sistivity. In the Harman method this same phenomenon is 
turned into a virtue. Consider a sample of a good thermo- 
electric with cross-sectional area A and length L suspended 
by thin wires in a vacuum of order 10° torr and through which 
passes a direct current I for long enough time so that a steady- 
state is reached. Assume no heat loss along the connecting 
leads nor any loss by radiation to the surroundings. The po- 
tential difference across the sample V, is made up of a purely 
ohmic voltage V, augmented by a thermoelectric voltage SAT 
which has its origin in the temperature gradient resulting from 
the Peltier effect (remember, the Peltier heat is generated at 
one end and absorbed at the other end and this sets up a 
thermal gradient across the sample). Using the definition of 
the figure of merit, z = S*/px, and some straightforward ma- 
nipulation, it is easy to show that the dimensionless figure of 
merit can be written as 





Hence, in principle, by measuring the adiabatic and isother- 
mal resistivities (voltage V_ with a direct current and volt- 
age V, with a chopped dc of equal magnitude) one can readily 
determine the dimensionless figure of merit ZT. 

In practice it is not possible to provide total thermal 
isolation of a sample from its surroundings and small heat 
losses along the wires supporting the sample and losses due 
to radiation from the end contacts and from the sample sur- 
face are inevitable. Correction terms which need to be ap- 
plied are given in a book by Goldsmid”*. With a judicious 
choice of the connecting leads and appropriate sample di- 
mensions the figure of merit at room temperature can be 
measured with an accuracy better than 2%. 


5. 


The Hall effect 


The Hall effect represents a powerful experimental tool 
which provides insight into the nature of transport, the den- 
sity of carriers, the type of conduction, and the carrier dy- 
namics. Measurements of the Hall effect have become not 
only a standard weapon in the armory of a solid state physi- 
cist seeking to understand transport behavior of solids but 
also an integral part of the quality control process in essen- 
tially every branch of solid state electronics. Although the 
thermoelectric figure of merit does not explicitly depend on 
the value of the Hall constant, the information gained from 
the Hall measurements is truly indispensable for a thorough 
understanding of the carrier transport in thermoelectric ma- 
terials. It is for this reason that I include here a brief descrip- 
tion of the Hall measurements. 

The Hall effect is a manifestation of the Lorentz force 
acting on a free charge carrier moving in a solid. As the car- 
riers are deflected to one edge of the sample (the top edge in 
the case of electrons and the configuration of Fig. 6a) they 
set-up a transverse electric field which opposes further de- 
flection and pile up of the charge. At a steady-state the trans- 
verse electric field, called the Hall field, exactly balances 
the Lorentz force. In experiments one normally does not 
measure electric fields but rather voltages. Writing the cur- 
rent density as j. = I/wt where I is the current and wt the 
cross-section of the sample, and applying magnetic field 
along the z-axis, the Hall voltage V,, becomes 


(9) 


Please note that it is the thickness t (the dimension in the 

direction of the magnetic field) and not the width w (the spac- 

ing between the Hall voltage probes) that enters in the Eq. 9. 
The quantity R defined as 


(10) 


is the Hall constant. Eqs. 9 and 10 are the essence of the 
Hall effect measurement. The direction of the transverse elec- 
tric field and therefore the sign of the Hall effect depend on 
the sign of the charge carrier. The established convention is 
that when q is negative (electrons), the Hall coefficient R is 
taken as negative. Fig. 6a, in fact, depicts the situation where 
the charge carriers are electrons. Furthermore, Eq. 10 facili- 
tates a straightforward determination of the carrier density. 
It holds very well for metals and strongly degenerate semi- 
conductors because here the carriers participating in the con- 
duction process occupy a narrow energy range on the order 
of «,T around the Fermi level. Having nearly the same en- 
ergy and therefore equal velocities, all carriers will be de- 
flected in a given magnetic field by approximately the same 
amount. In nondegenerate semiconductors, on the other hand, 
the carriers do not have the same velocities and they will be 
deflected in a magnetic field by different amounts. However, 
even in this case, the correction to Eq. 10 is small and can be 
accounted for by introducing a constant A on the right-hand 
side of Eq. 10 which never departs very far from unity. For 
instance, in a covalent solid where the carriers scatter on 
acoustic vibrations, the mean-free path is independent of 
energy and the constant A equals to 37/8 ~1.17. In ionic semi- 
conductors, at temperatures above the Debye temperature, 
the carriers scatter predominantly on optical phonons, and 
the constant has a value A = 1.11. 

By combining the Hall measurements with that of the 
electrical resistivity, one arrives at an important transport 


parameter known as the Hall mobility, ,, = o|R|. As dis- 
cussed elsewhere in this journal issue, the mobility is one of 





Figure 6a. 
Hall effect for a single band of free electrons. 


Figure 6b. 
Hall effect for a solid with a band of electrons and a band of 
holes. 
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the critical parameters which serves as a criterion whether a 
given material is a good thermoelectric or not. Large values 
of 4, are desired, and materials possessing a large value of 
the mobility are also expected to show significant changes 
in their transport behavior when subjected to an external 
magnetic field. 

If more than one group of carriers participate in the 
conduction process we have more than one Lorentz force 
and we must treat the steady-state as arising from the equal- 
ity of currents (magnetic field-generated deflection current 
balanced by the current due to the Hall field) rather than the 
equality of the Lorentz force and the Hall field. For a fre- 
quently arising situation of the mixed conduction (one band 
of electrons with 6,=nq, and R,=-1/nq, and a band of holes 
with O,=pqu, and R,=1/pq) a general expression valid for 
all field strengths is”’ 


wu (nu? — pu?) + u2u>(n— p)B 


2 
q (nu, + pu,) +u2u2(n-p) B 





(11) 


Obviously, when the material is perfectly compensated, i.e., 
when n=p the field dependence disappears. If, in addition to 
being compensated, the mobilities of the electrons and holes 
are the same, then the Hall coefficient becomes zero. In the 
limit of large fields and large mobilities, or more precisely, 
when 1 B>>1 and 1 B>>1, Eq. 11 reduces to 


l 


a 


lim Be 


(12) 


and the high-field Hall coefficient is a measure of the carrier 
imbalance. It is noteworthy that Eq.12 applies regardless of 
the shape of the Fermi surface or the nature of scattering. 

In principle, the configuration of Fig.6 coupled with 
Eq.9 or, in the case of two-band conduction Eq. 11, provides 
for a straightforward determination of the Hall constant. 
However, it should be remembered that the Hall voltage sig- 
nal is rather small, typically several orders below the ordi- 
nary resistive voltage generated in the same sample with the 
same current, and one must minimize systematic errors. 

The most critical issue in the Hall measurements is a 
probe alignment. Although one can align Hall probes on thin 
film structures very accurately using lithographic techniques, 
in bulk samples one usually is not so fortunate and probe 
misalignment is always present. In addition to the Hall volt- 
age, the misaligned probes inevitably pick up resistive volt- 
age which is proportional to the degree of misalignment. To 
eliminate the effect of probe misalignment one reverses the 
magnetic field and measures the Hall voltage again. The true 
Hall voltage is then obtained by subtracting the two read- 
ings, 
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ne V,,(B) - V,,(—B) 
2 





V, (13) 


In fact, it is a good practice to reverse not just the magnetic 
field but also a direction of the current and calculate the Hall 
voltage from the mean of all four readings. 

Another important consideration in Hall effect mea- 
surements is the length-to-width ratio L/w of the sample. In 
short samples with large width w (and current contacts over 
the entire width w) the measured Hall voliage would clearly 
be smaller than the value based on Eq. 9 because the Hall 
field is effectively shorted. It has been show” that in order 
to eliminate the shorting effect of the current contacts, the 
length-to-width ratio should not be smaller than 3. It should 
be realized that one does not solve the problem by merely 
making the current contacts smaller. It is true that the short- 
ing of the Hall voltage is then not as serious as before but 
this comes at a cost of making the current distribution highly 
inhomogeneous, perhaps even a greater problem. 

Equations describing the Hall effect, just as Eq. 1 for 
the electrical resistivity, assume that the measurements are 
made under isothermal conditions., i.e., in the absence of 
any thermal gradients. We have already pointed out that this 
is an unrealistic scenario when investigating electrical resis- 
tivity of thermoelectric materials with dc currents. The use 
of a chopped dc current, the best remedy identified then, is 
also the best approach in measurements of the Hall effect. 

For completeness, we note that the Hall effect can also 
be measured using the van der Pauw arrangement discussed 
in section 1b. Referring to Fig. 2, van der Pauw derived‘ the 
formula for the Hall mobility 


AR i413 
p 


6 
Uy a (14) 


where 6 is the thickness of the plane lamella, B is the mag- 
netic field, p is the resistivity, and AR.,, |, is the change of the 
resistance R,, ,, due to the magnetic field. Again, the clover- 
shaped sample geometry has a distinct advantage in provid- 
ing a relatively large Hall signal at low level of dissipation. 


Conclusion 


An increased sensitivity to environmental concerns 
with regard to waste heat and the deleterious nature of chemi- 
cal refrigerants and compressed gases has once again brought 
the development of new thermoelectric materials and devices 
to the forefront. Whether the search for new and more effi- 
cient materials will be successful hinges directly on our ability 
to measure accurately the fundamental transport properties 
of these systems. This is all the more true when it is realized 
that many of the materials being investigated today are from 





classes of solids whose transport properties are virtually un- 
known and potentially unusual and exotic in nature. Whether 
the observed effects are representative of the actual proper- 
ties of the material or are an artifact of an ill-fated experi- 
mental technique is a critical issue which every researcher 
must consider on a case-by-case basis. 

In this review I have attempted to outline the most 
popular and useful techniques for measuring the transport 
properties of thermoelectric materials, and have highlighted 
the advantages and drawbacks of each. As I have pointed 
out, the choice of a particular technique and its successful 
application depends on the type, size, and shape of the mate- 
rial, the temperature range of interest, the desired accuracy 
of the data, and, to a large extent, the particular researcher’s 
taste and experience. As new and more interesting materials 
are discovered and developed, their accurate characteriza- 
tion remains a challenge requiring serious attention, and it is 
hoped that this review will assist researchers in choosing and 
implementing the techniques best suited for their particular 
needs and tasks. 
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Research Notes 


Thermoelectric Power 
Studies at the Naval 
Research Laboratory 


by Alexander Ehrlich 


In response to the needs and opportunities represented 
by thermoelectrics, the Naval Research Laboratory has de- 
veloped a research thrust which is integrated in the Navy’s 
overall thermoelectric research program. This thrust includes 
the theoretical work of Singh et. al. and Reinecke et. al. which 
has been reviewed elsewhere in this volume. The first in- 
volves the calculation of electronic and phononic structure 
which are a central issue in the thermoelectrically important 
parameters making up Z. The latter work explores the po- 
tential for figure of merit enhancement by quantum confine- 
ment of the electrons in the thermoelectric material, i.e. us- 
ing very thin films or fine wires as the active thermoelectric 
material. 

The Naval Research Laboratory has developed a flex- 
ible state of the art experimental capability for measuring 
the various electrical transport phenomena needed to evalu- 
ate new materials for their potential use in thermoelectric 
applications. This has put NRL in a position to evaluate new 
classes of materials which appear promising, as well as be- 
come a center for evaluating materials developed specifically 
for these applications by a variety of institutions. This work 
is now being carried out by Drs. Wendy W. Fuller-Mora and 
Alexander C. Ehrlich in the Transport and Electronic Struc- 
ture Section, headed by Dr. Ehrlich and Dr. Valerie Brown- 
ing. A dedicated thermoelectric apparatus capable of mea- 
suring the Seebeck coefficient (thermoelectric power) and 
electrical conductivity from room temperature to four de- 
grees Kelvin, with data acquisition taking place in a 
semiautomated way as the temperature is swept between the 
limits of interest, has been developed. A second apparatus 
designed to obtain accurate measurements of thermal con- 
ductivity on much larger samples of a specific geometry and 
also capable of measuring S and o on these samples has been 
constructed and employed to measure metallic systems. It is 
currently being modified to achieve accuracy on the semi- 
conductors which constitute most of the practical thermo- 
electric materials. This apparatus is capable of measurements 
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at room temperature and above and also a computer collects 
data in a semiautomated manner. Dr. Kevin Stokes, an NRL 
post-doc, is also developing an experimental capability to 
measure thermal conductivity by a method known as the 3 
omega technique. This recently developed approach to the 
problem largely avoids the errors caused by unknown heat 
loss present, to a degree, in other approaches to the measure- 
ment of thermal conductivity. 

The existence of these apparatus’ has resulted in col- 
laborations with a number of laboratories including aca- 
demic, industrial and governmental institutions. For example, 
a number of examples of both filled and unfilled skutterudites 
compounds prepared by Slack et. al. at RPI have been char- 
acterized. With Marlow Industries of Dallas, Texas the pos- 
sibility of accurate, direct measurements of all the param- 
eters in Z, including thermal conductivity, has been explored. 
Dr. David Demske of the Naval Surface Weapons Center, 
White Oak is developing techniques for making quantum 
wires to achieve enhanced Z’s. The development is being 
coordinated with NRL to facilitate the process as well as to 
assure that the materials are geometrically and otherwise 
suitable for measurement as the development continues. 

An interdivisional effort to develop aerogels as ther- 
moelectric materials will be undertaken at NRL. Dr. Karen 
Swider, Dr. Debra Rolison (Chemistry Division), and Dr. 
Celia Merzbacher (Optical Sciences Division) will develop 
electrically conducting aerogels and characterize them by a 
variety of techniques. The goal is to employ the unique struc- 
ture of aerogels to enhance thermoelectric performance (e.g. 
by exploiting the lower dimensional structure of these ma- 
terials which may give rise to quantum confinement). Addi- 
tionally, a collaborative effort between Jerry Meyer (Optical 
Sciences Division) and Northwestern University research- 
ers Ketterson and Wong is also underway. This effort seeks 
to exploit quantum-confined carriers in bismuth for next- 
generation thermoelectric cooling devices. 





Thermoelectric 
Cooling Systems 


by Gordon Duval 


Thermoelectric (TE) cooling systems are currently less 
efficient than conventional vapor compression systems, and 
therefore, are not widely used. As part of the Chloro-fluoro- 
carbons Replacement program, Naval Sea Systems Com- 
mand (NAVSEA) has funded research into improving ther- 
moelectric material efficiency. Recent breakthroughs in this 
research could make TE attractive for future ships. System 
level studies have shown that the total ship electric power 
requirement for a decentralized TE cooling system with im- 
proved TE units, approaches the power required by a con- 
ventional centralized vapor compression system. Decentrali- 
zation also results in cost savings due to reduced ventilation 
ductwork, the elimination of an isolated fan room, and modu- 
lar construction. 

Thermoelectric air conditioning units (TEAC) were 
developed and installed on the USS DOLPHIN, AGSS 555 
in the early 1970s. Westinghouse Model 21E and 22K/L, 
one-half ton units (10 modules per cabinet) were used in the 
forward electronic space, the pump room and the engine 
room. Eventually module replacement became necessary. 
Having learned where the problem areas were, redesign was 
in order. Advanced design TEAC modules (Figure 1) have 
been developed and manufactured at Naval Surface Warfare 
Center Carderock Division (NSWCCD), Annapolis; and two 


prototype modules are currently aboard the DOLPHIN in the 
forward compartment (with four original modules). These 
modules utilize standard TE materials. They are 50% lighter 
and projected to be cheaper to manufacture than the original 
modules. An additional four modules will be manufactured 
and installed later in FY97, continuing the TE operational 
presence begun in 1971. 

The TE material research program was in the past di- 
vided into three separate tasks, improving the standard TE 
material, bismuth telleride; improving TE material configu- 
rations (both short range); and developing a new family of 
TE materials based on multilayer superlattices (long range). 
Current work at Massachusetts Institute of Technology (MIT) 
and Lincoln Labs has the potential to substantially increase 
the figure of merit (ZT). This effort has focused on the use 
of two-dimensional quantum well multilayer superlattices 
(QWSL) material systems. Theoretical models of these sys- 
tems have demonstrated a potential for large increases by 
using quantum well structures (atomic level layering of TE 
material) which constrain electron and phonon interaction to 
two dimensions. Lincoln Labs is working with MIT to fabri- 
cate these multilayer systems through molecular beam epi- 
taxial (MBE) techniques. Experimental measurements in 
FY96 have confirmed ZT values greater than 2.3 for these 
PbTe quantum wells. It should be noted that for bulk PbTe 
material, ZT=.6 (approx.) is the best measured. Work will 
continue in FY97 to reduce barrier thickness, improve ther- 
mal conductivity measurements, and obtain higher three di- 
mensional ZT values. 





Figure 1 


Water-to-air thermoelectric cooling module. NSWCCD prototype. 
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